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Abstract 

These lecture notes give an introduction to a number of ideas and methods that 
have been useful in the study of complex systems ranging from spin glasses to D-branes 
on Calabi-Yau manifolds. Topics include the replica formalism, Parisi's solution of the 
Sherrington-Kirkpatrick model, overlap order parameters, supersymmetric quantum me- 
chanics, D-brane landscapes and their black hole duals. 
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1 Introduction 



1.1 The case for complexity 

If you are a string theorist, chances are that you think of a complex structure as something 
squaring to minus one, rather than as something consisting of many intricately interacting 
degrees of freedom. These lecture notes, however, are concerned with the latter. 

There is a deeply rooted belief in the natural sciences that the more fundamental a theory 
gets, the less important complexity becomes. In particular there has traditionally been an 
almost unquestioned assumption that physics at subatomic scales must be intrinsically simple, 
and that complexity is only relevant in the macroscopic or living world. Certainly there is 
a decrease in complexity when going from cells to proteins and from there to amino acids, 
atoms and finally elementary particles. It would seem logical that any step further down this 
reductive tree — for example explaining the origin of elementary particles and their properties, 
accounting for the entropy of black holes or understanding the initial conditions of the universe 
- should involve only structures and concepts of the utmost simplicity and elegance. 

By now considerable theoretical and experimental evidence has accumulated suggesting 
this is wrong. In particular, whenever quantum mechanics and gravity pair up — as in cos- 
mological dynamics on the largest space and time scales, the determination of effective low 
energy parameters in theories with compact extra dimensions, and the stringy microphysics of 
black holes — complexity appears not only unavoidable, but seems to some extent essential. 
On the theoretical side I am referring here more specifically to the highly complex, fractal-like 
iterated structures arising in eternal inflation, the perplexing complexity of the string theory 
landscape, and the closely related complexity of wrapped D-brane systems that, among many 
other applications, have given us the first quantitative explanation of the entropy of black 
holes. 1 On the experimental side, the most influential development has been the gathering of 
abundant cosmological precision data: detailed cosmic microwave spectroscopy providing con- 
vincing evidence for the slow roll inflation hypothesis, and the measurement of an impossibly 
tiny yet nonzero vacuum energy density, of just the right magnitude to allow self-reproducing 
resonances to crawl out of gravitationally collapsed dust at precisely the right time to see the 
onset of cosmic acceleration. Whether we like it or not, these developments, and the absence 
of plausible alternative explanations for these and other fine-tuning conundrums, have added 

Except for very specific items, references will be given at the end of this section and in subsequent sections. 
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significant credibility to the idea that some fine-tunings may have no other explanation than 
being tautological post-selection effects in a vast, eternally inflating multiverse, scanning a 
huge ensemble of low energy effective field theories, only a tiny fraction of which will allow 
structures to evolve of sufficient complexity to be capable of organizing summer schools. This 
idea requires no fundamentally new concepts of any kind, works by exactly the same semi- 
classical mechanisms as ordinary inflation, eliminates some of the worst fine-tuning problems 
in physics, provides in principle a precise, quantitative definition of naturalness, and appears 
to fit naturally into a UV complete, unifying framework incorporating all known principles of 
physics, that has no external parameters yet spawns all the necessary complexity, out of simple 
equations following unambiguously from an even simpler idea: the quantization of strings. 

The price to pay for this remarkable feat is that it forces us to consider structures and con- 
cepts far more diverse and complex than those it seeks to explain, without immediate prospects 
of directly observing them. The Standard Models of elementary particles and cosmology is 
arguably a lot simpler than any of the string compactifications that could produce it; for ex- 
ample it could be F-theory compactified on the Calabi-Yau fourfold hypersurface in weighted 
projective space CP ' ' 8 ' ' 806 , which boasts about three hundred thousand deformation 
moduli, up to seventy-five thousand D3-branes, and a flux-induced superpotential specified 
by a choice of almost two million integers multiplying an equal amount of independent period 
integrals [1, 2]. According to simple estimates, this compactification, assuming all moduli can 
be stabilized, gives rise to more than 10 1QJ intricately interconnected flux compactifications. 
Similarly, a macroscopic extremal Reissner-Nordstrom black hole looks considerably simpler 
than the corresponding weakly coupled D-brane systems one is led to consider in microscopic 
computations of the black hole entropy. For example a wrapped D4 producing a modest 
extremal black hole of say the mass of the sun easily involves D4-branes with 10 100 moduli 
subject to superpotentials specified by 10 100 flux quanta. Finally, even in simple toy model 
landscapes, the large scale structure generated by eternal inflation is infinitely more complex 
than what we observe in our universe. 

Although the importance and urgency of a better understanding of these matters is clear 
to many, there has been widespread reluctance to face this kind of stringy complexity directly. 
Instead, the dominant approach has been to try to isolate particular phenomena of interest, for 
example by studying simplified compactification models like local Calabi-Yau spaces, effective 
field theory models, non-disordered toy models of the landscape, D-brane models for black 
holes restricted to charge regimes where the dominant contributions to the entropy are highly 
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structured by symmetry (Cardy regime), or coarse grained asymptotic parameter distributions 
in ensembles of compactifications allowing similar complexity-minimizing limits. There are of 
course obvious excellent reasons to follow this reductionist approach. However, there are also 
excellent reasons to try to probe complexity itself: 

1. In contrast to supersymmetric AdS compactifications, a specific choice of compactifica- 
tion data does in general not correspond to a superselection sector in cosmologies with 
a positive vacuum energy, as quantum and thermal fluctuations will force changes in ge- 
ometric moduli, fluxes and topology. Thus, any complete, nonperturbative description 
of a sector of string theory rich enough to describe our own universe should be able to 
encode not just one compactification, but the full space of internal space configurations 
that can be dynamically reached in one way or another. 2 Since various topological tran- 
sitions between for example Calabi-Yau manifolds are known to be perfectly sensible 
physical processes in string theory, the space of such interconnected configurations is 
likely to be huge. If we imagine for a moment that there exists a complete holographic 
description of eternally inflating cosmologies in string theory, say in the form of a field 
theory living at future infinity, then this field theory would somehow have to encode 
the full complexity of those googols of geometries — it would literally be a theory of 
everything, and it couldn't be anything less. 

2. There is often striking organization, universality and elegance emerging in disordered 
systems, as has become clear over the past decades in studies of spin and structural 
glasses, neural networks and other complex systems. Conventional notions of symmetry 
are largely irrelevant in such systems, but other, equally powerful structures and their 
associated order parameters appear in their place. This includes hierarchical cluster 
organization of the state space, replica symmetry breaking and overlap order parameters. 
Uncovering these led to highly nontrivial exact solutions of various models of complex, 
disordered systems. Much like symmetries in ordered systems, these structures also 
determine to a large extent the dynamics and other physical properties. They have a 
wide range of applicability in fields as diverse as condensed matter physics, neuroscience, 
biology and computer science, and have led to practical applications such as new efficient 
algorithms for optimization, data mining and artificial intelligence. Thus, rather than 
an annoyance to be avoided, complexity can be the essence, and the key feature to focus 

2 It is not necessary that such configurations also support metastable "vacua" . 
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on. The lessons learned from these investigations are bound to be useful in the context 
of complex systems in string theory and cosmology. 

3. Conversely, string dualities such as holography may provide a new and useful approach 
to the general theory of complex systems. What makes this particularly promising com- 
pared to other instances of applied holography is that the natural large N limit is also 
the limit one is a priori interested in here. Concrete examples are given by supersym- 
metric branes wrapped on compact cycles in Calabi-Yau manifolds, which exhibit many 
of the typical characteristics of mean field models of glasses. 

1.2 Contents and goals of these notes 

In these lecture notes, I will give an introduction to some of the concepts and techniques which 
have been useful in studies of complex systems of many degrees of freedom, with intricate, 
disordered interactions. Although I will discuss ideas and techniques developed in the theory 
of glasses as well as ideas and techniques developed in string theory, there will be little or no 
discussion of applications of one to the other. For this I refer to the work that will appear in 
[3, 4, 5], and in which any new idea that might be present in these notes originated. The focus 
will be on the basics, meaning material usually assumed to be known in most of the recent 
specialized literature. Rather than to give a comprehensive review, I will treat a number of 
specific topics in a more or less self-contained and hopefully pedagogical way, to avoid variants 
of step 7 of [6] , and to allow the reader to learn how to actually compute a number of things 
rather than to just get a flavor of the ideas. This comes at a price of having to leave out many 
interesting and important topics (including some that were discussed to a certain extent in my 
actual lectures at TASI). In particular, despite its importance in the motivation given above, 
I will say almost nothing directly here about the landscape of string compactifications, and 
focus instead on space localized D-brane systems arising as the weak coupling description of 
charged black holes. This has three additional reasons: (i) I have already written extensive 
lecture notes on string vacua [7], and wanted to avoid overlap altogether, (ii) The internal 
space geometry of the space-localized D-brane systems describing black holes is identical to 
that of space-filling D-brane systems describing compactifications, and the techniques used to 
analyze the former can directly be transported to the latter. The advantages of looking first 
at space-localized branes are furthermore numerous: They are conceptually and technically 
easier, allow arbitrarily large charges and hence a proper thermodynamic limit, and have 
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holographic dual descriptions as black holes, providing effective "experimental measurements" 
of thermodynamic quantities like the entropy. (Hi) It prevents, for the time being, potential 
diffusion of confusion from partially unresolved conceptual problems of quantum gravity in a 
cosmo logical setting (as reviewed by Tom Banks at this school [8]) into what should be an 
exposition of well-understood and diversely applicable methods. 
More specifically, I will cover the following topics: 

1. An introduction to the theory of spin glasses, in particular Parisi's solution of the 
Sherrington-Kirkpatrick model, using the replica formalism. This was the first non- 
trivial energy landscape to be studied and understood in detail in physics. Special 
attention is given to Parisi's overlap order parameters, which allow to detect a non- 
trivial equilibrium state space structure without having to be able to explicitly know 
those equilibrium states. A generalization of this order parameter for arbitrary quan- 
tum systems is proposed at the end of the section. Some other complex systems such as 
the Hopfield model for memory and learning are briefly discussed. Other, non-replica 
approaches such as the cavity method and Langevin dynamics are important but not 
treated in any detail in these notes. 

2. An introduction to supersymmetric quantum mechanics, where simple but powerful 
concepts such as the Witten index make it possible to compute exact quantum ground 
state degeneracies of highly complex systems. The computation of nonperturbative 
lifting effects due to landscape barrier tunneling as well as relations to Morse theory 
are discussed in some detail. Finally, the map between nonsupersymmetric Langevin 
dynamics and supersymmetric quantum mechanics is introduced, and its appearance in 
glass theory is outlined. 

3. An introduction to the low energy quantum mechanics description of D-branes wrapped 
on compact cycles in Calabi-Yau manifolds. These are the prime examples of complex 
systems in string theory, and excellent models to introduce techniques ubiquitous in 
the field of string compactification. Again I start from scratch, assuming only some 
background knowledge in elementary differential geometry and the basics of D-branes 
in string theory. I describe in detail the low energy reduction of the D4-brane wrapped 
on a high degree 4-cycle. This results in a supersymmetric quantum mechanics of high 
complexity, which nevertheless is manageable thanks to the underlying supersymmetric 
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geometry. If the D4 charge is of order N, the brane allows of order N 3 worldvolume mag- 
netic fluxes, inducing an extremely complicated superpotential for its order N 3 moduli, 
and leading to an enormous energy landscape with exponentially many minima. Never- 
theless it is possible to explicitly construct vast numbers of exact supersymmetric critical 
points, by giving the critical point condition the interpretation of "capturing" holomor- 
phic 2-cycles. At large DO-charge, this method becomes inadequate, but at the same time 
it becomes easier to give good estimates for the number of ground states. A formula for 
the index is derived by directly applying the machinery of Lagrangian supersymmetric 
quantum mechanics, treating the worldvolume fluxes as momenta canonically conjugate 
to angular variables. This reproduces the well known continuum approximation at lead- 
ing order, but also provides all corrections in terms of differential euler characteristics 
of critical point loci. Finally, bound states with mobile DO-branes are introduced and 
counted in a rudimentary fashion, a microscopic formula for the D4-D0 is derived, and 
some open ends and generalizations are summed up to conclude. 

4. An overview of the construction of the zoo of multicentered supersymmetric black hole 
bound state solutions that are the holographic duals of the D-brane energy landscapes 
discussed in the previous section. 

1.3 Some general references 

For the topics covered in these lectures, I will give references along the way. For the topics 
mentioned above but which will not be treated in any detail in what follows, here are a few 
general referencesthat may be useful: [9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 
24, 25, 26, 27, 28, 29, 30, 22, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 7, 41, 42]. 

2 Spin glasses and beyond 
2.1 Introduction 

A spin glass [43, 44, 45] is a system of localized spins with disordered interactions. An exam- 
ple is copper sprinkled with manganese atoms at random positions. The sign of the spin-spin 
interaction potential between the manganese atoms oscillates as a function of distance and is 
relatively long range, so we get effectively randomized mixed ferromagnetic and antiferromag- 
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netic interactions. Some spin pairs will want to align, while others will wish to anti-align. If 
among three given spins, two pairs want to align and one pair wants to anti-align, or if all 
three pairs want to anti-align, one pair will not get what it wants and the triangle is called 
frustrated. The presence of many frustrated triangles typically leads to exponentially many 
local minima of the energy, and thus to a landscape. 

This landscape is highly complex; the problem of finding the ground state in even the 
simplest spin glass models is effectively intractable. Nevertheless, as we shall see, some models 
of spin glasses have exactly solvable thermodynamics, thanks to a remarkable underlying 
mathematical structure which greatly simplifies the analysis but nevertheless leads to a very 
rich phase space structure. Its most striking consequence is ultrametricity of the state space. 
An ultrametric space is a metric space with a distance function d that satisfies something 
stronger than the standard triangle inequality, namely 



What this means is that all triangles are isosceles, with the unequal side the shortest of the 
three. Ultrametric spaces appear in various branches of mathematics; for example the p-adic 
distance 3 between rational numbers is ultrametric. For nice introductions to the mathematics 
of ultrametric spaces, see chapter 9 of [46], and [47]. Ultrametric spaces are also familiar in 
biology: If we define the distance between two current species as how far back in time one 
has to go to find a common ancestor, this distance is trivially ultrametric. Indeed, picking a 
turkey, a monkey and a donkey for example, the distance turkey-donkey and turkey-monkey 
is the same, while the distance monkey-donkey is the shorter one of the three. Alternatively, 
one can define a species distance as the degree of difference of DNA. Since this will be roughly 
proportional to how far back a common ancestor is found, this distance will again be essentially 
ultrametric. The same is true for any representative set of proteins. This is used with great 
success to reconstruct evolutionary trees (for a recent account see e.g. [48]). The lesson here 
for us is that ultrametricity is equivalent to the possibility of clustering objects and organizing 
them in a hierarchical tree, similar to the evolution tree of species. 

In the case of spin glasses, the points in the metric space are the different equilibrium states, 
that is the different ergodic components or superselection sectors in which the standard Gibbs 
measure breaks up in the glass phase. We will describe this in detail in the subsequent sections, 

3 The p-adic norm of a rational number is \q\ p = p~ k when q can be written as q = p k m/n where m and n 
contain no powers of p. The p-adic distance is then d p (a, b) = \a — b\ p . 
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Figure 1: An ultrametric tree. Distances between the points in the top layer are set by the vertical distance 
to the first common ancestor node. 

here we give a qualitative picture. The distance between different equilibrium states is given 
by some measure of microscopic dissimilarity, for example the sum of local magnetization 
differences squared. Remarkably, in mean-field models such as the Sherrington-Kirkpatrick 
model [49], this space turns out to be ultrametric [50, 51, 52]: In the thermodynamic limit, 
with probability 1, triples of magnetization distances are isosceles. Equilibrium states get 
hierarchically organized in clusters. The cluster division is independent of the distance measure 
used, as long as it is statistically representative, meaning zero distance implies equal states. 
Figure 2 shows the results of a simulation illustrating this phenomenon. 

What is remarkable is that such a nontrivial tree structure appears in a purely static 
setting, without any apparent underlying evolution process creating it. 4 Instead, the role of 
evolutionary time is played by temperature or more generally energy scale: At high tempera- 
tures there is only one ergodic component of phase space. When the temperature is lowered 
below the spin glass phase transition, the phase space starts breaking up into distinct parts 
that no longer talk to each other, separated by effectively unsurpassable free energy barriers. 
Initially these valleys in the free energy landscape will all still look very similar, but upon 
further cooling, their mutual magnetization distance grows, and eventually they in turn will 
start breaking up in components separated by free energy barriers. This branching process 
proceeds all the way down to zero temperature, generating an ultrametric tree of states. Any 
pair of ground states defines a "time" on this tree, proportional to their ultrametric distance. 

4 To be precise, the remarkable thing is the existence of a nontrivial ultrametric structure, with multiple 
branchings. Trivial ultrametricity, meaning all triangles are equilateral, is very easy to realize in infinite 
dimensional spaces. For example all pairs of randomly chosen points on an infinite dimensional sphere are 
with probability f at the same distance from each other. 
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Figure 2: [5] Dendrogram plots and overlap matrices for the SK model with N = 800 spins. We used parallel 
tempering Monte Carlo [53, 54] to reach thermal equilibrium, with 50 replicas at equally spaced temperatures 
between T = 0.1 and T = 1.2. Based on overlaps of 100 configurations sampled with separation of 100 sweeps, 
clustered using Mathematica. Results are shown for T/T c = 1.2,0.86,0.55,0.12. Red = maximal positive 
overlap q = +1 (i.e. minimal distance), white = overlap q = 0, dark blue = maximal negative overlap q = — 1. 
Although the data is too limited to draw firm conclusions, the plots are suggestive of hierarchical clustering. 

As we will see, ultrametricity plays an important role in the exact solution of mean field 
spin glass models, and determines many of their static and dynamic properties. It plays in 
this sense a role similar to that of symmetries and symmetry breaking in physics. In fact, its 
emergence is closely related to the breaking of an auxiliary symmetry, called replica symmetry, 
which naturally arises in the description of disordered systems. The order parameter (or rather 
the order function) capturing the symmetry breaking pattern exhibits ultrametricity as one 
of its characteristic features. We will discuss this extensively in what follows. 

For more realistic, local (short range) spin glass models in three dimensions, for which the 
mean field approximation is poor, there is no consensus on the physical presence or relevance of 
ultrametricity, and different schools exist with different favorite models and results, the most 
widespread being the droplet model of [55]. On the other hand, in systems where little or no 
locality is present, such as neural networks or combinatorial optimization problems, mean field 
methods and ultrametricity are ubiquitous [43]. A recent numerical study of ultrametricity in 
simple Ising type models of different degrees of locality can be found in [56]. 

For a review of ultrametricity in physics, see [57]. General criteria for the appearance 
ultrametricity were formulated in [58]. Introductory texts on spin glass theory include [43, 
59, 60, 61] 
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In the following we will make the above qualitative description more precise, introduce 
the new order parameters and explain the mysterious but powerful replica trick and the 
phenomenon of replica symmetry breaking associated with ultrametricity. 



2.2 The Parisi order parameter 
2.2.1 Spin glass models 

The most basic models for spin glasses are formally similar to the classical Ising model, with 
N spins Si = ±1 interacting according to a Hamiltonian of the form 5 



In the case of the Ising model, we have Jjj = — J for nearest neighbors (ij) and zero otherwise. 
At high temperatures the magnetization M = 4 rrii is zero (wij is the average value of Sj). 
In dimension 2 or higher, below a critical temperature T c , the magnetization acquires a nonzero 
value and the Z 2 spin flip symmetry is spontaneously broken. In the mean field approximation, 
T c = J and M is a solution of the mean field consistency equation M = tanh(M J/T). 

In 1975, to model spin glasses, Edwards and Anderson [62] proposed to study the above 
Hamiltonian, still with nearest neighbor couplings J^, but now drawn randomly and indepen- 
dently from a Gaussian distribution with zero mean and standard deviation J. Although at 
finite N everything in this model will depend on the actual values of Jy, in the thermody- 
namic limit N — > oo, intensive quantities such as the (spatially averaged) magnetization or 
free energy density are self-averaging: they become equal to their J^-averaged value with unit 
probability. This makes computations possible. Due to the disorder, the magnetization M 
vanishes at all temperatures in this model. Nevertheless, there exists a critical temperature 
below which the spins freeze, in the sense that they acquire locally preferred directions. This 
disordered frozen phase is the spin glass phase. The magnetization is clearly not a good order 
parameter to detect this situation. Instead Edwards and Anderson introduced a new order 
parameter: 



In the statistical mechanics literature, the common convention is to put a minus sign in front of the sum 
over the spins, so ferromagnetic couplings are positive. We will not do this in these lecture notes. 




(2.2) 




(2.3) 
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Here the TOj are the local magnetizations in the equilibrium state the system finds itself in, 
which in an actual physical setting can be thought of as the time averaged value of the spin s;. 
For the Ising model mi is independent of i, so q E A = M 2 . The Edwards- Anderson parameter 
is related to the experimentally measured magnetic susceptibility: x = (1 — q E j^)/T. 

Though the meaning of time averaging is physically clear, we have not specified any dy- 
namics for this model, so we need a purely static definition for m^. This is subtle. It is not the 
canonical ensemble average (sj) obtained from the Gibbs probability measure p(s) oc e~^ H ^ s \ 
as this average is trivially zero due to the Z 2 symmetry of H. The same issue arises already 
for the Ising model, but there it is clear what is going on: below T c the Gibbs probability 
measure splits in two "superselection" sectors, characterized by opposite values of the mag- 
netization. The two sectors can be separated by switching on a small background magnetic 
field h, which shifts H — >• H + /^X^ 5 *' lifting the Z 2 degeneracy and eliminating one of the 
two sectors, depending on the sign of h. We can then define the magnetization for each sector 

as mfct = lim^oi lirnjv->.oo(si)fc- 

For spin glasses this does not work; since M = 0, switching on a constant h will not 
lift the degeneracy. Switching on an inhomogeneous field hi tailored to the m; of the frozen 
equilibrium state we want to single out would of course do the job, but it is impossible to 
know in advance what the required profile is going to be; it could be anything, and it will 
be different for different values of Jy. To make things worse, there may be many distinct 
equilibrium states, all with different values of the and again there is no way of telling 
without actually solving the system. There is no standard symmetry and no order apparent 
in the system at any temperature. What "order" parameter could possibly distinguish between 
those phases? 

These questions were considerably sharpened and then answered through the study of an 
even more simplified spin glass model, proposed by Sherrington and Kirkpatrick [49] in 1975. 
The Hamiltonian is as above except that now the couplings are nonzero for all pairs of 
spins, not just nearest neighbors; they are all independent Gaussian random variables with 
zero mean and variance J 2 /2N, 6 

JijJkl = ^jySikSjl ■ (2.4) 

Overlines will denote averages over the throughout these notes. The model is completely 

6 The variance is scaled with N such that the typical size of the interaction potential at each lattice site 
remains finite in the thermodynamic limit: (J^. Jys J ) 2 = J 2 /2. The factor of 2 is added for consistency with 
standard conventions, in which usually Jij = Jji is imposed, which we will not do. 
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nonlocal; spatial geometry has no meaning in this setup. It is exactly solvable, but the solu- 
tion turned out much more interesting than Sherrington and Kirkpatrick originally imagined. 
The simple solution they originally proposed was correct at high temperatures, but mani- 
festly wrong near T = 0, where it predicted for instance negative entropy densities. After 
considerable effort, the correct solution was finally found in 1979, in a seminal breakthrough 
by Parisi [50, 51]. (It took another 25 years before the solution was rigorously proven to be 
correct [63, 64].) In the course of the process, Parisi uncovered the proper order parameters 
to fully describe the spin glass phase, as well as a hidden "statistical" symmetry group whose 
breaking they parametrize. We will describe the solution in section 2.3, but first we describe 
the physical meaning of the order parameter. 7 To do this we must introduce the notion of 
thermodynamic pure states. 

2.2.2 Pure states 

The probability measure characterizing the Gibbs state is 

PoW = |e-^«, Z = Xy/>*«. (2.5) 

s 

For a given probability measure p, we denote the expectation value of an observable A(s) 
by (A) p . For the Gibbs state we usually just write (A) G = (A). We say p satisfies cluster 
decomposition if correlation functions of local observables factorize in the thermodynamic 
limit for almost all points. That is, for any finite r, in the limit iV — > oo, 

\Ai x Bi 2 ■ ■ ■ Ci r ) p = (Ai 1 ) p (Bi 2 ) p ■ ■ ■ {Ci r ) p + Ri x i 2 ...i r , (2-6) 

where the remainder R is negligible on average: liniAr^oc -^j- ^i r ..i l-^n -vl = 0- Here A iy Bi, 
. . . , Cj are local observables, like for example A^ = Sj or Bi = Sj_iSjSj + 3. In particular the 
clustering property implies that intensive quantities like the free energy density have definite, 
non-fluctuating values. In local theories (2.6) is equivalent to the property that correlation 
functions factorize in the limit of infinite spatial separation. 

The Gibbs measure for the Ising model satisfies cluster decomposition above the critical 
temperature, but not below: at high temperatures we have (siSj) = = (si)(sj) when — > 
oo, but at low temperatures we have instead (siSj) = M 2 ^ (si)(sj) =0-0. However we can 

historically it was first obtained as a formal mathematical object in the replica formalism. The physical 
interpretation was given later in [65]. 
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canonically decompose the Gibbs measure into a direct sum of probability measures which 
do satisfy cluster decomposition. Thus, for the Ising example, we have below the critical 
temperature 

p G (s) = lp + (s) + ip_(a) , (2.7) 

where p±(s) = lim^-to limjv->oo ~e~^ H ^ s '' Jrh ^ Si ' . Sets of finite measure according to p + have 
zero measure according to p_ and vice versa. The measures p± do satisfy cluster decomposi- 
tion: (siSj) + = M 2 = (si) + (sj) + . The superselection sectors described by p± are called "pure 
states" in the statistical mechanics literature. 

Rigorously defining pure states for general systems is subtle. For a discussion in the 
context of spin glasses we refer to page 89 (appendix 1) of [66]. What is known rigorously is 
that any probability measure pc defined for a system of infinite size, which locally behaves 
like the Gibbs measure (i.e. it gives relative probabilities proportional to e~^^ H for finite size 
fluctuations), can always be uniquely decomposed into pure states p a as 

Pg(s) = ^2w a p a (s) . (2.8) 

Here w a > is the Gibbs probability of being in the pure state a. The pure states satisfy the 
clustering property and cannot be further decomposed. This gives them an implicit definition. 8 

2.2.3 Overlap distributions 

As mentioned earlier, for disordered systems, it is in general impossible to explicitly find 
the actual decomposition into pure states. But granting the decomposition exists, one can 
formally define the overlap q aj 3 between pure states as 

Qafi = ^2(Si) a (Si)fi ■ (2.9) 

i 

As a special case, notice that q aa is nothing but the Edwards-Andersen order parameter q^A 
defined in (2.3). Although not at all obvious, it turns out [67] that q aa is independent of 
the pure state a (and of the disorder realization J^), so q^A is actually an invariant of the 

8 More intuitively, in physical systems, pure states correspond to distinct ergodic components, individually 
invariant under time evolution and not further decomposable into smaller time invariant components. In 
ergodic components the ergodic theorem implies that the time average of any observable equals its ensemble 
(phase space) average. 
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system, depending only on T/J. 9 Using this, we see that the overlap is closely related to the 
Euclidean distance between states: 

d(a, (3) 2 = — ^2 - m^) 2 = 2(q EA - q a/3 ) . (2.10) 

i 

At this point the overlap matrix may still seem like an abstract, incomputable quantity. 
However, consider the overlap probability distribution: 

P{q) = ^2waWp6(q- q a p) . (2.11) 

a/3 

This is the probability of finding an overlap q a p = q when one samples the Gibbs state. The 
wonderful thing is that this quantity is actually computable without any knowledge of the 
actual decomposition into pure states. It can be rewritten purely in terms of the Gibbs state 
as the overlap distribution for two identical replicas of the system, with spins and s^: 

P(Q) = (s(q-^E l sf ) sl 2) )) n=2 , (2.12) 

where n = 2 means we are considering two replicas: 

<^ ( V 2) )>_ 2 = ^ ^e-^ (1) )-^ (2) )A( S «, S ( 2 )). (2.13) 

s 

To prove this, one shows the moments are equal. Consider for example its second moment: 
(q 2 ) = JdqP{q)q 2 

aft i j 

ap ij 

- Hb ™)'L- 

9 Since the SK model is completely nonlocal and does not a priori distinguish any site j relative to a given 
site i, we can think of the permutation symmetry acting on the site indices i as a gauge symmetry, similar 
to the diffeomorphism group in gravity (where Jjj is the analog of the spacetime metric). Then the local 
magnetization distribution V{m) = A S(m — m ia ) = 5(m — m ia ) encodes all gauge invariant information 
based on the "vevs" rrii a . It is shown in [67] that again this quantity is independent of the state a. So we can 
say that there is no gauge invariant distinction between different pure states based on just the magnetizations; 
all equilibrium states look the same as far as non-fluctuating quantities are concerned. 
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Figure 3: From [68]. Disorder-averaged overlap distribution P{q) for the SK model at various temperatures, 
from T — 0.95 on right to T = 0.30 on left in steps of 0.05 (T c — 1). The dotted lines represent delta-functions, 
localized at q aa = Qea- 

To get to the third line we used the clustering property (2.6), and for the fourth line we used 
(2.8) in reverse to express everything in terms of Gibbs state expectation values. The higher 
moments are treated analogously. 

Equation (2.12) in principle allows us to compute or approximate P(q) by standard meth- 
ods, for example by using Monte Carlo simulations sampling the Boltzmann-Gibbs distribu- 
tion. It captures the presence and properties of pure states very well and is therefore a good 
order parameter for the spin glass phase. If only one pure state is present, namely the Gibbs 
state itself, as is the case at high temperatures, then q u = (Jea = due to the Z 2 symme- 
try, and P(q) = S(q). If the system freezes and splits into one Z 2 pair of pure states, with 
local magnetizations rrii and —rrii, then for both states the self-overlap is qEA — ^Yli m l 
and P(q) = \5{q — qEA) + \&(q + 9ea)- If there are many pure states, there will be many 
delta-functions, and if there is an infinite number of them, P(q) becomes continuous. 

Computing P(q) analytically for a given coupling matrix is not possible. It is however 
possible to compute its average P(q), by a variety of methods that can also be used to compute 
more basic thermodynamic quantities; the most prominent ones are the replica, cavity (TAP), 
Langevin dynamics and supersymmetry methods [43, 59, 69, 70, 71, 72]. In these lecture 
notes we will focus on the replica method. To already get an idea of where all this is heading, 
a result of such computations is shown in fig. 3, displaying the disorder-averaged P(q) for 
various temperatures below T c . The function is smooth except for a delta-function peak at 
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the highest value of q, which is the self-overlap q aa = qEA- Its presence is a consequence of 
the state and disorder independence of qEA, leading to a term {^2 a w^)6(q — qEA) in (2.11). 
The nonvanishing of Yla w a i n turn indicates the presence of states of finite weight w a . The 
continuity of P(q) shows that there are in general multiple pure states, but it does not imply 
Pj(q) is continuous at fixed J; in fact because some w a are finite, Pj{q) must have multiple 
delta peaks, but with J-dependent locations, so integrating over J smooths them out. Indeed, 
unlike the magnetization or the free energy, Pj(q) is not a self-averaging quantity; it fluctuates 
between different realizations of the disorder. 

2.2.4 Ultrametricity 

One can also define more refined overlap distributions, such as the overlap triangle distribution 
P{qi,Q2,q3) = ^2w a w f }W-yS(q 1 -qf} 1 )S(q 2 -q 7a )5(q3-q a i3), (2.14) 

a/3'y 

which similar to (2.12) can be written in terms of the overlap distribution of three independent 
replicas in the Gibbs state. Restricting to triples with positive overlaps, 10 one finds [52] 

P{QuQ2,qs) = 2 J rf </W) p (<?i)<Kgi - £2) %2 - qs) 

+ \ ( P ^ P ^ ®( qi ~ q ^ ^ ~ 93 ) + P ermut ations^ , (2.15) 

where 9 is the step function and P(q) is as in (2.11). This manifestly exhibits ultrametric- 
ity, with the first term encoding equilateral and the last three terms more general isosceles 
triangles. 

2.3 Replica solution of the SK model 
2.3.1 The replica trick 

We now turn to the actual computations in the SK model, first using the replica formalism 
[62, 49, 50, 51]. The Hamiltonian is, as stated earlier: 

H jI s ) = ^ JjjSjSj , (2.16) 

10 This is necessary because without this restriction, the Z2 symmetry trivially destroys the ultrametric 
structure: A sign flip of the spins in state a will flip the sign of two out of three overlaps of an (a/37) triangle, 
and clearly this does not preserve isoscelesness. This also explains the structure in the blue regions of fig. 2. 
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where i,j = l,...,N,8i = ±1 and Jij is drawn randomly out of a Gaussian distribution with 
mean and variance i.e. 

p{Jij) oc e'^i . (2.17) 

Other distributions would give equivalent results; only the first two moments matter. It is 
impossible to solve this system exactly for arbitrary given couplings — just finding the 
absolute minimum energy configuration is already an NP-hard problem [73]. Fortunately 
the usual intensive thermodynamic quantities of interest are self- averaging, that is they are 
independent of the random matrix with probability going to 1 in the thermodynamic limit. 
Thus we can compute them by computing averages. For example the free energy density is 

F = Fj= fUdJ^piJ^Fj, Fj = -±-\ og Zj , Zj= e-^'H. ( 2 -!8) 

J ij sd{±l} N 



Crucially, the average log Zj has to be computed after taking the logarithm. Such an aver- 
age is called quenched, the disorder represented by the is called quenched disorder, etc. 
Computing the average first, i.e. on the partition function itself, is called annealed averaging. 
Physically, this corresponds to a situation in which the couplings themselves are fluctuating 
variables. The annealed average is much easier to compute than the quenched average, but it 
is not the situation we are in; in a real spin glass for example the couplings are determined by 
the positions of the impurity atoms in the host crystal, which vary randomly over space but 
do not fluctuate in time, or if they do, on much longer time scales than the spins fluctuate. 

To deal with the quenched average, one can use the replica trick. This is based on the 
observation that 

Z n _ i Q 

logZ = lim = Z n \ n=0 , (2.19) 

n-s>o n on 

together with the fact that for positive integers n, the average of Z n is just an annealed average 
of the disorder coupled to n replicas of the original system, and therefore easy to compute. 
Of course, the subtle part is the "analytic continuation" to n = 0, which is hard to make 
rigorous. 11 However, its effective power in a large range of applications is undeniable, so the 



11 A perhaps somewhat more precise version starts from the observation that Z n is an entire function of n, 
with Taylor expansion Z n = e nl °z z = J2T=o jr(\ogZ) k valid for all n, so = 1 +n\ogZ+ ■ ■ ■ for all n. The 
idea is then that the term linear in n can be extracted by considering arbitrary integers n and computing Z n 
as an expansion in n. However this is in principle ambiguous: for example 1 + sin(7m) is entire and evaluates 
to 1 for all integers, yet it has a nonvanishing first order term in its Taylor expansion. For a deeper analysis 
see e.g. [74, 75]. 
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lack of a rigorous framework is no doubt more a failure of our arsenal of rigorous frameworks 
than of the idea itself. 

Introducing replica indices a,b = 1, . . . , n, the coupling-averaged n-fold replicated partition 
function is, putting J — 1, 

W = j dJij e~ NJ l J2 e-P J » s « si . (2.20) 

se{±i} nN 

Repeated indices are summed over and we absorb the Gaussian normalization factors in the 
measure dJij. At this point, the replicas do not interact with each other. However, because 
they all couple to the same Jij, integrating out the disorder induces effective interactions 
between the replicas: 

W = e^'«"k*"*. (2.21) 

s£{±l}" N 

The positive sign in the exponent means the interaction is attractive. Now, reverse-mimicking 
how replicas got coupled by integrating out the site coupling matrix J^, we decouple the 
lattice sites by integrating in a "replica coupling matrix" Q a t,: 




(2.22) 



The S'-sum here runs over the n replica copies only; lattice indices no longer occur. The sum 
over the lattice sites has simply produced an overall power of N. This is awesome because 
it means that in the thermodynamic limit iV — > oo, we can evaluate the integral over Q in a 
saddle point approximation. It all boils down now to finding the critical points of 

HQ) = 4 l°g Z{Q)> Z(Q)= E e^ QahSaSh - (2-23) 

Se{±i}™ 

Denoting the dominant critical point(s) of J-{Q) by Q+, we thus get Z n = e -P N:F (Q+)_ jf we 
can find the saddle points for general n, we are done. Using the trick (2.19) and noting that 
consistency with Z° = 1 requires J r (Q*)| n =o — 0, we obtain the free energy density: 

F=^-F(Q*)\ n=0 . (2.24) 
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To summarize, what we have done is trade summing over the lattice for summing over replicas, 
which allows us to do a saddle point computation of the quenched average. Before proceeding 
to find the solution, we pause to ponder the meaning of the matrix Q ab that apparently 
captures the large N behavior of the model. 



2.3.2 Meaning of Q 



(lb 



So far the physical meaning of the matrix Q ab is obscure, although the conspicuous notational 
similarity with the overlap matrix q a p introduced in (2.9) suggests the two are related. We 
will now show that this is indeed the case, in that its structure in the limit n — > will produce 
for us the Parisi order parameter P(q) defined in (2.11), or more precisely P{q)- To see this, 
first observe that the saddle point equations ^'(Q*) = can be written as the self-consistency 
equation 



1 



Q,a b = (S a S b ) Q ^—- S a S b e^ Q * M . (2.25) 

Hence Q ab equals the overlap in replica space. Now consider the expression (2.12) for P(q): 

P(q)= E & %-^E s ^)- (2-26) 

si,s 2 e{±i} JV i=i 

This depends on the random couplings J^. We want to average it but are facing a problem 
similar to the problem we had when we wanted to compute the average of log Z: now the 
problem is the explicit appearance of J in the factors 1/Z. We deploy again the replica trick, 
this time based on 

1 N 

P ^ = iTo E e ^ Ea=1 * [Sa] %-]vE s i4)> (2-27) 

se | ±1 }nJV i=l 

where we have singled out the first two replicas to appear in the delta function. Indeed, the 
right hand side can for n > 2 be written as 

Z n-2 e~^ H ^ +H ^ S(q - j- E 44) , (2-28) 

si,s 2 G{±l} Ar »=1 

formally reducing to (2.26) when continued to n — > 0. Since in (2.27) there are no longer any 
denominators, we can compute the disorder average by simple Gaussian integration. 
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Using (2.27) and manipulations similar to those leading up to (2.22), one computes the 
finite moments (q k ) = J dqP(q)q k . Dropping l/N suppressed terms: 

/2 / 2 \ N— k / 2 \ k 

^5e{±i}" ' ^5e{±i}" ' 

f / _ p ^QabSaS b \ * 

= nm '«( £ 4wT**) • < 2 - 29 > 

^ v Se{±i}™ wy 7 

where we have absorbed in the measure dQ a normalization factor 1/ J dQ e~^ NJ7 ^\ which 
ensures (q°) = 1. The functions £(Q) and J-"(Q) were defined in (2.23). We may replace Si £2 
by S a Sb with a 7^ b and average over replicas, since there is no distinction between replicas. 
The resulting expression has the advantage of being the same at distinct saddle points related 
by the permutation symmetry of the replicas. In view of (2.25), the saddle point evaluation 
then gives 

[Q*] V 1 *#> [Q*] 6^1 

The sum over is a sum over distinct replica permutation symmetry orbits, i.e. it is non- 
trivial only if there are different saddle points not related by the permutation symmetry of 
the replicas, as may be the case when there are additional symmetries beyond the replica 
permutation symmetry. For each such orbit Q*, Wq^ denotes its relative weight. In the case 
at hand, we do have an additional Z 2 spin flip symmetry in the system. If this is the only 
other source of saddle degeneracy, we have Wq+ — \. In general we conclude: 

W) = - lim W Q* E % " > ( 2 - 31 ) 

[Q*] m 

and if there is no degeneracy besides the Z2 one (as one would expect generically) , this becomes 

P{q) = - lim V \5{q - Q* lb ) + \5{q + Q* lb ) . (2.32) 

Ml 

Thus the pure state overlap distribution will be entirely determined by the saddle point 
solution Q* continued to n — > 0. To obtain the probability to find an overlap q, all we need to 
do is compute the fraction of entries Q+u that are equal to q. The weird looking minus sign 
is not a typo: it must be there because Ylb^i 1 = (ft — 1) — >■ — 1 when n — > 0. Still, counting 
the number of entries Qu, in a x matrix may cause a feeling of unease in the reader. In 
the following section we will try to make sense of this, and conclude that, paradoxically, the 
properly continued 0x0 matrix Q ab in fact has infinitely many degrees of freedom! 
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2.3.3 Replica symmetric solution 

Finding the critical points of F(Q) for general n is still a nontrivial task. Some obvious conse- 
quences of (2.25) are Q* a a — 1 and the fact that it is a positive definite symmetric matrix, but 
to make further progress, one must make an ansatz for the form of the solution. The simplest 
one, which was the one used by Sherrington and Kirkpatrick [49], is the replica symmetric 
(RS) ansatz, which is the unique ansatz leaving the permutation symmetry unbroken (we will 
drop the explicit * subscripts from here on): 

Qab = uSab + q(l ~ 5 a b) ■ (2.33) 
So Q a b = u if a = b and Q ao = q if a ^ b. The relation (2.32) then gives 



P(q') = -(6(q'-q)+6(q' + q)). (2.34) 

Thus the RS ansatz is equivalent to assuming not more than one Z2 pair of pure states. If 
q = 0, there is just one pure state and the system is in the paramagnetic phase. 

We wish to extremize F(Q) = \Qa\~~ \ 1°§ e^® abSaSb . From the ansatz (2.33) we obtain 

Ql b = nu 2 + n(n - l)q 2 , e ^Q^s h = e ^(n( u - g)+g (E a &) 2 ) _ 

Assuming without loss of generality q > 0, we linearize the term involving S a in the exponent 
with the transformation 

Se{±i} n V n< l J 5e{±i}" 

= —L= fdze'i [2cosh(/3z)] n (2.37) 
V 27i q J 

= 1 + n^= / ^e"^log[2cosh(/3z)] + 0{n 2 ) , (2.38) 
\/27rg J 

whence 

d n F\n=o = j{u 2 - q 2 ) - |(« -q)- i-^= J dze~£ log [2 cosh^)] . (2.39) 
This is to be extremized with respect to u and q. Extremizing u is trivial: 

u = 1 , (2.40) 
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reproducing the result Q aa = 1 we arrived at earlier directly from (2.25). We set u = 1 in 
what follows. There is no simple closed form solution for the saddle point value of q, so let us 
consider limiting cases. In the high temperature limit (3 — > 0, we expand in powers of (3: 

d n T\ n=0 = -{ (1 + q 2 ) ~ ^ + 0(/3 2 ) • (2.41) 

The extremum is at q = 0: As expected for high temperatures the system is in its paramagnetic 
phase. The free energy and entropy densities at q = are, from (2.24): 

r - l0g2 S = f3%F = hg2-^. (2.42) 



(3 4 ' ~ - ~o- 4 

The log 2 corresponds the two-fold spin degeneracy at each lattice site. 

Below a critical temperature T c , q = ceases to be the thermodynamically stable saddle 
point. To find this transition point, we expand (2.39) for small q 

a j~i log2 t f3((3 2 -l) 2 (3 5 3 4 

d n JF\ n= Q = + q - y q 6 + 0(q ) . (2.43) 

The coefficient of q 2 changes sign when (3 = 1, and a new saddle point with q w /3 — 1 > 
becomes the thermodynamically stable one. 12 This signals the spin glass phase transition. 
Recalling we chose units such that the coupling variance parameter J = 1, we conclude 

T C = J. (2.44) 

This seems all fine, and is consistent with Monte Carlo simulations. However, in the low 
temperature limit (3 — > oo something awkward happens. We have 




d n F\ n=0 = - v / 7 (3 - \ — + 0(/3- 2 ) . (2.45) 



7T 



The relevant saddle point isatg = l — y-^. This means the self-overlap q^A of the pure 
states approaches 1, as expected. But the free energy and entropy density are, to leading 



12 What "stable" means in the limit n — > is actually rather subtle [76]. The reader may notice for example 
that in the case at hand, the stable saddle points are local maxima in g-space, whereas ordinarily in physical 
parameter spaces, stable equilibria are those that minimize the free energy. This "inverted" rule finds its 
origin in the formally negative dimension of the fluctuation modes of Q a b in the limit n = 0. For example, 
the expression U = — ^ a b Q 2 ab is manifestly positive definite for positive integers n. It has a positive definite 
Hessian (the n 2 x n 2 unit matrix) for all n and the Q = extremum of U is stable. Nevertheless, inserting the 
RS ansatz (2.33) gives U(q) = 1 — q 2 . It is in general nontrivial to do a full stability analysis in the replica 
formalism. For more discussion see [76, 591. 
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Figure 4: Tree and matrix representation of a Parisi matrix Q a b- Different colors and different tree connection 
heights correspond to different values of Q a b, with darker red representing larger values. In this example, 
n = 90, K = 3, {m ,TO 1 ,m 2 ,m3,m 4 } = {1,3,15,45,90}, {",(70,91,92,93,94} = {1,0.9,0.6,0.4,0.1,0}. 

order 

- + -0.798 , 5=-^- w -0.159. (2.46) 

7T ZTC/j Z7T 

The free energy disagrees with numerical simulations, which indicate F ~ —0.76 at T = 0. 
More dramatically, the negative entropy clearly does not make any sense. Therefore the 
solution must be wrong at low temperatures. The replica symmetric ansatz (2.33) is apparently 
too restrictive. The replica symmetry must be broken. 

2.3.4 Replica symmetry breaking and Parisi matrices 

We want to relax the replica symmetric ansatz (2.33), which was Q ab = 5 a b + q(l — 8 a b). 13 

The effective replica Hamiltonian H R = — ^ ^2 a f,{jf Si s a s l) 2 appearing in (2.21) de- 
scribes an attractive interaction between the replicas: It is energetically favorable if the repli- 
cas line up. Competing with this is the fact that lining up means less phase space. Below T c , 
the energy gain is more important than the entropy loss and we get a nonzero overlap expec- 
tation value between the replicas. In the RS ansatz the degree of correlation is assumed to be 
the same for all pairs of replicas. Relaxing this assumption means considering the possibility 
that groups of replicas form clusters that have larger overlaps amongst themselves than with 
other clusters. The maximally symmetric situation in such a scenario corresponds to the case 



13 We put the diagonal part of Q a b equal to 8 a b here and it what follows. This is justified by observing that 
the diagonal part of Q a i, decouples from the off-diagonal part in (2.25), and that for any saddle point Q aa = 1. 
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in which the clusters are indistinguishable from each other, in particular equal in size and 
mutual overlaps. In this case there is also no absolute distinction between the replicas, only 
relative to some fixed other replica. 

To find the thermodynamically relevant critical points of F(Q) (defined in equation (2.23)), 
it is therefore natural to consider [50, 51] a minimal permutation symmetry breaking scheme 
S n —> S n j m x S m , by splitting up 14 the n replicas in clusters of size m, and assuming the 
overlap between distinct replicas within one cluster to be qo and the overlap between replicas 
in different clusters to be q\ < qo. Refining the replica labeling by where ai = 

1, . . . , n/m labels the clusters and a = 1, . . . , m the replicas inside each cluster, this translates 
to 

QaiacMo = ^aifei^a 6o + QO ^a^^aobo + <?1 ^aibi > e ab = 1 _ $ab (2.47) 

(so t a b = 1 if cl 7^ b and if a = b). The idea is then to substitute this new ansatz in J-(Q), 
extract the 0(n) term, and look for saddle points by varying q , qi and m. 

The clustering process can be iterated, by breaking up the clusters into smaller clusters 
and then those in turn into even smaller ones. When there are K distinct nontrivial cluster 
sizes, this is referred to as replica symmetry breaking at level K (i^-RSB). 

To write things out explicitly, put m = 1 and denote the size of the smallest nontrivial 
cluster by mi, the next larger one by m 2 and so on, up to mx+i = n. Let be the overlap 
between replicas within a cluster of size m^+i (excluding those contained in an even smaller 
cluster). Labeling replicas by a = clkclk-i • • • ai^o where a« = 1, . . . , labels the clusters 

771% 

of size jrij, the ansatz can be written as 

K 

Qa K —a ,b K —b = $ab + q% $a K b K ' ' ' Sg i+1 b i+1 £ ai bi (2.48) 

i=0 

= A i 5 *Kb K ■ ■ ■ <W>i , Ai = -qt > 0, (2.49) 

i 

where qi = 1 if i < and qi = if % > K. For example for K = 2 this becomes 

Q a2CL\aQ,b2b\bo "(12^2 ^d\bl ^ao^O *?0 "0262 ®ai fei ^(M)bo Q 1 "0262^0161 Q2 ^0262 

= (1 - 90)^262^01 61 <W + (<?0 ~ ?l)£a 2 &2<Wi + (?1 - <?2)<W 2 + <?2 • 

14 When m does not divide n, it is of course not possible to split the replicas in equal blocks of size m. The 
idea is to consider values of n and m for which it is possible, and then to analytically continue the resulting 
expressions, treating m, like n, as a continuous parameter. In this section we only consider the integer n case. 
The continuation will be discussed in section 2.3.5. 
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Reading this from right to left corresponds to zooming in to finer structures of the overlap 
matrix. An example with K = 3 is shown in fig. 4. 

The hierarchical block structure is equivalent to having an ultrametric structure in the 
overlap matrix. In other words, it can be organized as a tree, as illustrated in fig. 4. Defining 
an ultrametric reference distance r a b between two replicas to be, say, the size m of the smallest 
cluster to which they both belong, we can reformulate the Parisi ansatz simply as the statement 
that the overlap only depends on the distance: 

Qab = q{r ab ) , (2.50) 

where q(r) is an arbitrary function (encoding the qi). 

We want to evaluate the replica free energy J-(Q) = \Q 2 ah — ^ log e~ QabSaSb defined in 
(2.23). The first term is straightforward: 



nh \ i — n / 



K 

(2.51) 



ab i=0 

The interpretation of the coefficient (m i+1 — m^) multiplying each qf is clear: it is the number 
of replicas having overlap with some fixed reference replica. The second term in J-(Q), 

/3 2 a b 

proportional to log Z(Q), where Z{Q) = e^® abSas , is more interesting. Using (2.49), we 
have 



Z = ^exp y^Ai ( S "K~a ) 



2-[ 

(2.52) 



The squares in the exponential can be linearized by Gaussian transforms similar to (2.36). 
The resulting expression initially involves Gaussian integrals over many variables z~a K ... ail but 
they can be evaluated iteratively in clusters starting from the smallest cluster and integrating 
up to larger and larger distance scales (i.e. larger % / larger / smaller q^). The result of 
this little exercise is conveniently and suggestively expressed in terms of convolutions with the 
Green's function of the heat equation, 

G ««^ exp [-Cl< (2 ' 53) 
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Figure 5: The cluster size function m(q) for the K = 3 example of fig. 4, and the corresponding tree. The 
arrow indicates the direction of the RG flow. 

producing Z as the outcome of the following recursion: 

Z (z) = ^V^ s = 2cosh(/3z), (2.54) 

s 

Z i+ i(z) = dz' ' G {qi „ qi _ l{ (z - z') , k i=^, (2.55) 

Z = j dzG qK {z)Z K+1 {z). (2.56) 

The powers hi arise from identical copies of Gaussian integrals, the copies corresponding to 
different values of the sub-index a t = 1, . . . , The first step in the recursion is easy enough: 

/ „2 \ mi 

Zi(z) = e- (1 - 9o) (2cosh(/3z)) , but past this point one has to resort to expansions or nu- 
merical evaluation. As a simple check of the above result, notice that when we remove the dis- 
tinction between the clusters, i.e. qx = • • ■ = <h — Qo (but keeping g_i = 1, qx+i = 0), we have 

/ B 2 , \ m K+l 

Gfa-^iz - z') = 5{z - z') for 1 < i < K, so Z = J dzG qo (z) ( e V<^)(2coBh(0*))) 
Recalling m K+ i = n, we see this correctly reproduces the replica symmetric formulae of section 
2.3.3. 

The recursion (2.55) can be thought of as an exact Wilsonian renormalization group action, 
evolving from the UV (small clusters) to the IR (large clusters), rescaling the number of degrees 
of freedom in jumps by factors ki = m i+1 /mj. We can make the equations more familiar 
looking by defining a function Z(q, z) for all q G [0, 1] such that Z{qi — e, z) = Z i+ i(z). The 
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equations (2.55) and the identification of G as the Green's function of the heat equation show 
that we can take this function Z(q, z) to be the single spin sum 2 cos\i((3z) at q — 1, and then 
evolve it down in q according to the heat equation d q Z = — ^d 2 Z, with a jump Z — > Z kl 
whenever q crosses a gj. As the latter transformation is nonlinear, it is more convenient to 
work with g(q,z) = log Z(q, z), since on g it acts linearly: g — > kig. (The price to pay for 
this is that the heat equation becomes nonlinear.) To implement this without reference to the 
indices i, we define a function m(q) (illustrated in fig. 5) as 

m(q) = size of clusters with overlap > q . (2-57) 

So in particular m(gj) = and d q m(q) = ^(mj — m i+1 ) £(g — gj). Theng(g,,2) = \ogZ(q, z) 
is the solution to 



v q9 = -\ 



dg\ d 2 g 
dz I dz 2 



d 1 dm 
dq m dq 



V q = — --—, (2.58) 



evolving down from q = 1 to q = with initial condition 



and final identification 



lim g(q, z) = log (2 cosh(/3z)) , (2.59) 



\ogZ = \img(q,0) . (2.60) 

q\0 



The "gauge connection" A q = (d q m)/m keeps track of the increase in the number of degrees of 
freedom, while the nonlinear equation in g is equivalent to the heat equation. The connection 
may be gauged away by redefining / = ^-g, so d q f = -^Dqg, for which 



df = _1 
dq 2 



m{q) \dz + dz 



dfY , d 2 f 



(2.61) 



We thus see that the iterative system is nothing but a particularly simple example of the exact 
renormalization group equations, with q the analog of the scale and f(q,z) the analog of the 
Wilsonian action. 

For some purposes it is more convenient to take the scale variable to be some arbitrary 
auxiliary parameter A, parametrizing the (q,m) staircase (fig. 5) by continuous functions q(X) 
and m(A). Letting A evolve from A = to A = 1, with 

(q, m) | A=0 = (1, 1) , (q, m) A=1 = (0, n) , (2.62) 
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and denoting A-derivatives by a dot, the RG equation becomes 

'df\ 2 , d 2 f 



f = - q - 
J 2 



m(A) 



dz J dz 2 



/(z)| A=0 = log(2cosh(^)). (2.63) 



Obviously this introduces a reparametrization gauge symmetry. Putting everything together, 
we get 

ijfy, m] = ^1 + jf dm q 2 ^j - ±f(z = 0) | A=1 . (2.64) 

The remaining tasks are (i) make sense of the analytic continuation to n = and (ii) extremize 
the functional J-[q,m]. 

2.3.5 n ->• with RSB 

Given that 1 = m < mi < m 2 < - - - < m^ < rriR+i = ^ in the Parisi matrix construction, 
it is not obvious how to properly continue to n = < 1, to say the least. To get an idea of 
what could constitute a sensible continuation, consider the overlap distribution (2.32), which 
in i^-RSB becomes, taking into account the crucial overall sign of (2. 32), 15 

K 

P(q) = lim y (mi — m i+ i) 5(q — qi) = lim d q m(q) . (2.65) 

8=0 

For this to make sense as a probability density we need 1 = m > mi > m 2 > • • • > mx+i = 
n = 0, i.e. the above cluster size inequalities must be inverted, and the obviously can 
no longer be integers; we will allow them to be arbitrary real numbers. Equivalently, m(q) 
must become an increasing function from the unit interval to the unit interval, to make d q m 
a proper probability density. This also inverts the interpretation of m(q) from (2.57) to 

limm(g) = probability of finding an overlap < q. (2.66) 

71— >0 

In general the function m(q) may have continuously increasing parts, which can be thought 
of as the K — > oo limit of the discretized construction. This is illustrated in fig. 6. Recalling 
the discussion of section 2.2.3, getting a smooth function should indeed not surprise us: even 
if individual realizations of the disorder produce a discrete overlap structure like the one 
depicted on the right of the figure, if this structure itself is sensitive to the disorder, averaging 
over the will smooth out the discreteness. 



15 To avoid having to duplicate everything, we restrict again to positive overlaps here, i.e. we condition the 
probabilities on q > 0. 
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Figure 6: (Fictitious) example of (m, g)-line after continuation to n = 0. The arrow again indicates the 
direction of the RG flow. The diagram on the right is a representation of the ultrametric state space tree 
for some choice of the disorder J^-. The vertical lines are separating the branches of the tree. Different 
compartments at the top (the leaves) correspond to different pure states a, and their width represents their 
probability weight w a . The value of q where two branches separate sets their mutual overlap (so the taller the 
wall, the smaller the overlap). Thus, going up along the q-axis, m{q) increases in discrete jumps. But different 
choices of lead to different barrier structures, and averaging produces the smooth overlap distribution on 
the left. (The residual step up is due to the disorder independence of the self-overlap qEA-) 

This suggests what we should aim for. Looking back at the equations (2.62)-(2.64), we 
see that in fact, the required structure is very naturally obtained by simply taking n — > 
in (2.62) while keeping (2.63) and (2.64) the same. Indeed, m(A) becomes now a decreasing 
function while q(X) remains decreasing. Everything else stays the same. In particular the 
intrinsic ultrametric structure is preserved. Notice that now the quadratic terms J* A=Q dm q 2 
in (2.64) become negative definite, with q = being the maximum. This is not an indication 
of thermodynamic instability, but rather the converse, as mentioned already in footnote 12. 
In the limit n — >■ 0, thermodynamically stable states correspond to maxima of the free energy 
J^iQ) instead of minima. 

Admittedly, the above is not a particularly solid justification for the proposed continuation 
to n — 0. One can probably do better, but I will not try to do this here. 16 

16 The continuation procedure is somewhat reminiscent of continuation in the p-adic numbers: A p-adic 
integral over a region with p-adic norm between 1 and n is a finite sum with p n terms, and it naturally 
continues to a p-adic integral over a region with p-adic norm between n — and 1, which is an infinite 
series. Similarities with p-adic and adelic cumbers were explored in chapter 9 of [46] and [77, 78]. However to 
reproduce the above one still needs to do a formal continuation p — > 1. 
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2.3.6 The Parisi solution 



To summarize, the Parisi solution 17 of the SK model is (going back to the gauge A = 1 — q), 
in all its glory: 

"\=o 



max 

m(q) 



1 



dm q' 



q=0 







(2.67) 



where m(q) is an non-decreasing function on the unit interval satisfying the boundary condi- 
tions 

m(0)=0, m(l) = l, (2.68) 
and f m (q, z) is the solution to the following flow equation running from q = 1 to q = 0: 



df 
dq 



m(q) 



df 

dz 



+ 



dz 2 



9=1 



log(2cosh(/32)) . 



(2.69) 



Although this may still look somewhat unwieldy, in practice a simple trial family of m(q) 
functions (for example stepwise constant functions with just a few steps K) can be treated 
easily and already lead to excellent approximations of the exact result. For example [51], 
for K = 1 one gets for the ground state free energy F = —0.765 and entropy S = —0.01, 
already significantly better than the RS result (2.46). For K = 2 this improves further to 
F = —0.764, S = —0.004. In particular we see that the entropy converges rapidly to the 
exact S = 0. Detailed computations at various temperatures can be found in [68]. A recent 
stability analysis of the low T solution is given in [79]. 



2.4 Beyond spin glasses 

This is as far as we will go in describing the explicit solution of the Sherrington-Kirkpatrick 
model. We have left many topics untouched, in particular the various other ways in which 
the model can be solved (for which we refer to [43, 59] as a starting point). We hope we have 
made it clear that the hierarchical organization of the state space was the central structure in 
the story. This appears to be typical for disordered mean field (i.e. nonlocal) models. Other 
popular examples of such models include the Ising and spherical p-spin models [80, 81, 61] . The 

17 To compare to [51]: x, h there equals m resp. z here, and the apparent extra terms there are generated by 
the initial flow from q = 1 to q = qEA during which m(q) — 1, which is included here but not there. This part of 
the flow is easy to solve because it is basically exactly the heat equation: f(q, z) = ^-(1 — q) +log(2 cosh(/3z)), 
as is easily checked. 
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latter has continuous variables Sj constrained to lie on an iV-dimensional sphere ^\ sf = N, 
with a Hamiltonian H = ^ ir ..j p Ji^-ipS^ ■ ■ ■ Si p , where the J are random couplings. It is 
known [61] that when p = 2, the replica symmetric solution is exact, while for p > 3, the 
1RSB solution is exact. In the limit p — > oo, the model becomes equivalent to the "random 
energy model" [80] and can be solved exactly and explicitly [81]. One can interpolate between 
spherical spin models and Ising spin models such as SK by switching on an additional non- 
random potential V = J2i ^( s l ~~ 1) 2 > where p = 2, A — > oo reproduces the SK model. Various 
generalizations to higher spins, SU(N) spins, rotors, quantum spins etc have been studied 
too, for example in [44, 82, 83, 84]. 

As mentioned in section 2.2.1, such mean field models are not expected to be a good 
approximations for physical, short range spin glasses in spaces with a moderate number of 
dimensions (including three), and it is not clear to what extent the hierarchical organization 
of the mean field models is a useful zeroth order expansion point to understand the physics of 
such systems. But there are other systems of interest, in physics and beyond, where nonlocal 
interactions occur naturally, and where the conceptual framework developed for the mean 
field theory of spin glasses has turned out to be very useful. This includes the theory of 
neural networks, either as models for the brain or for various forms of machine learning, like 
the way Google figures out all the time what you really mean. One of the oldest and most 
famous such models is the Hopfield model [85], which is effectively described by a Hamiltonian 
H = J^ij Jij s i s j, where the couplings Jy encode synaptic strengths and Sj = ±1 encodes 
the on/off firing state of the i-th neuron. Stable firing states of the neurons correspond to 
minima of H. iV-component patterns & with £j = ±1 are stored according to a learning rule 
Jij — > Jij — ^i^j (or variants thereof), which energetically favors the stability of neuron states 
Si = £j. After M patterns are learned, the couplings are thus = jjj — Y^!=i CfCj 4 - However, 
as is well known, learning does not equal remembering. Memory retrieval dynamics in this 
context is modeled by a standard Monte Carlo type relaxation process, which tries to minimize 
the energy. Remembering a learned pattern £ M then means that the change in coupling 
effectively has given rise to a stable minimum s* at or very near the pattern Clearly this 
model shares many characteristics with mean field spin glass models, and mean field spin glass 
theory has indeed been very useful in the analysis of these models. For example in [86] it was 
shown that with the above learning rule, the memory undergoes a phase transition when it 
tries to learn too much: for M < M c with M c rs IAN it remembers practically everything, 
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while for M > M c it gets completely confused and remembers practically nothing. 18 In this 
framework, the phase in which learned patterns are remembered is similar to the spin glass 
phase (many highly metastable states), hierarchical clustering in state space gets interpreted 
as categorization and organization of memories in concepts, and so on. A collection of some 
of the original papers along these lines can be found in [43] . 

Another important 19 branch of human activity where the ideas of spin glass theory have 
found fruitful applications is optimization algorithms for NP-hard problems. Canonical ex- 
amples of such problems are the satisfiability problem and the traveling salesman problem, 
which find incarnations in practical problems going from optimal chip design to the problem 
of finding string vacua satisfying certain conditions [40]. For an early review we refer again 
to [43]. Recent developments include for example the relation of clustering phase transitions 
to effective hardness of optimization problem instances and ways to exploit these insights in 
devising new efficient algorithms [89, 90]. 

2.5 Beyond spins 

The statistical mechanics community is fond of Ising spins, and for a good reason: there is a 
lot in the real world that can be naturally modeled by large collections of interacting binary 
degrees of freedom. String theorists on the other hand face the sad reality that the complex 
systems that are natural in their world are rarely as simple. On the other hand, Parisi's 
idea of using the overlap distribution to characterize the emergence of different pure states / 
ergodic components / superselection sectors in an intrinsic, quantitative and computable way 
is an attractive one. Thus, when pondering for example the state space structure of complex 
D-brane systems [5] or de Sitter space [4] , one is led to the question what the analog is of the 
pure state spin overlap q a p = 4 J2i( s i)a( s i) P defined in (2.9) for general quantum systems. For 
the more general systems of interest the degrees of freedom could be fluxes, scalars, fermions, 
KK modes, or whatever. Since typically there is no canonical scalar product between such 
degrees of freedom, a naive direct generalization of the overlap will not work. 

A simple construction that is fully quantum and makes sense in general, and that is 
equivalent to the standard spin overlap in the Ising spin case, is the following. Since we 
want to consider systems for which we can at least formally take a thermodynamic limit, 

18 It turns out that the key to retain the ability to remember new things is forgetting old things [87]. 
19 As measured by the 2 x 10 4 citations of e.g. [88]. 
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we assume we can randomly sample 20 N independent "representative" degrees of freedom Xi, 
i — 1, . . . , N of the system, with the thermodynamic limit of infinite system size corresponding 
to N — > oo. For example consider a scalar field theory living in a large box of volume V. 
We can subdivide the large box into smaller identical boxes of some fixed volume v, and then 
uniformly sample at some fixed sampling density a subset of N of the smaller boxes, which we 
label by i = 1, . . . , N. Finally, in each of the sampled boxes i, we compute the average value 
of the scalar field, and we call this Xi. Another example would be a U(M) matrix quantum 
mechanics describing a wrapped D-brane system. Here we may sample a uniform subset of iV 
of the M 2 matrix degrees of freedom. A final example is again an Ising spin glass, for which 
we sample N spins. 

The state of the complete system is assumed to be described by some density matrix p. 
Expectation values are computed as usual as (A) p = Tr(pA). For each degree of freedom i 
there is a corresponding reduced density matrix p % obtained by tracing out all other degrees 
of freedom, so 

Tr % (pVL t ) =Tr (pA t ) . (2.70) 

Here Hi is the Hilbert space associated to the degree of freedom x^ (for the scalar field example 
this would be Hi = £ 2 (M), whereas for (quantum) Ising spins it is Hi = C 2 ), and A+ is an 
observable referring to Hi alone (so for the scalar this could for example be Ai = x 2 , or 
Ai = —id x i). 

A natural definition of overlap between two states p a , p@ is then 

1 - 

i=\ 

When the system consists of Ising spins, this reduces to the standard magnetization (or spin) 
overlap, up to a rescaling and a shift. For indeed, the general density matrix for a spin 1/2 
degree of freedom can be written as p = | (1 + m • s), where m = (s) p is the magnetization, 

20 Sampling is not necessary for the construction, but it spares us the trouble of having to know and work 
with the complete set of degrees of freedom, which may lead to UV divergences or may not even be known. 
Moreover the overlap q a p and corresponding spin glass order parameter P{q) could equally well have been 
defined by sampling say every other spin instead of all spins; ultrametricity implies that the results would 
have been identical. In general this need not be true of course, so the notion of overlap may depend on the 
choice of degrees of freedom that are sampled. 
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whence Q a/ 3 = jj £\ |(1 + ' ^Jg)- 21 For continuous degrees of freedom this can also be 
written in terms of the associated Wigner densities on phase space: 



N 



= ^ J2 27lh I drdpW* a (r,p) Wj(r,p) , (2.72) 
i=i J 

where W^(r,p) = ^ J ds e ips / h (r — s/2\p % a \r + s/2). In the classical limit this becomes (after 
suitable coarse graining to kill off highly oscillatory modes) the overlap of ordinary phase 
space probability densities, and we can think of Q a/ 3 as the fraction of degrees of freedom 
(x l ,pi) found in the same elementary phase space cell (of size 2nh) when randomly sampling 
from two states a and /3. 

We can mimic (2.9) even more closely, by writing the reduced density matrix p l as the 
expectation value (E l ) of the operator E l defined by 

El s = Vl =r V Xi=s = \ Xi = r)( Xi = s\ , (2.73) 

where V Xi=r is the projection operator on the eigenspace Xi = r. Then we get 

1 - 

Q°f>=JfYlY l ( Ei rs)a(K)fi- (2-74) 
i=l rs 

Finally, for some states p, an analog of the clustering property (2.6) will hold, for others 
not. We define this analog to be the property that correlation functions of "local" observables 
factorize in the thermodynamic limit for almost all evaluation points. That is, for any finite r 



lim — ^ 



(AiiBi 2 ■ ■ ■ C ir ) p — {A i± ) p (B i2 ) p ■ ■ ■ (C ir ) f 



0. (2.75) 



By "local" we mean here observables A$ referring to Hi alone (we could relax this definition 
by allowing dependence on other, "nearby" degrees of freedom, but this will be good enough 
for our purposes). 

With these definitions set up, we can proceed in complete parallel with the Ising spin 
glass case, and define "pure states" p a , decompositions p = w a p a of non-pure states (such 
as p = e~@ H for glassy systems at low temperatures), and their overlap distributions P(q) = 
Yliap w a w ^{Q~Qai3)- Whenever p fails to satisfy cluster decomposition (signaling a nontrivial 



21 This is valid for quantum Ising spins. Classical Ising spins are trivially obtained from this by restricting 
to the up/down eigenstates of s z , yielding Q a p = A (1 + q a p) with q a p as defined in (2.9). 
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phase structure), P(q) becomes nontrivial. For example for an isotropic quantum ferromagnet, 
the low temperature pure states will be labeled by a vector rh such that for all i, p^ — 
|(1 + m • s), and the overlap distribution P(q) is uniform on the interval [ 1- 2 m2 , l+ ™ 2 ] . 

Moreover P(q) will again be computable in principle for any known state p, even if we do 
not know how to describe its decomposition into pure states. Explicitly, the analog of (2.12) 
that allows this is 

P{q) = (H<1-W J2irs ElMsr) ) n=2 , (2.76) 

where again the subscript £ n = 2' means we consider two replicas of the system, with density 
matrix p £g> p, and Ei depends on the first replica's degrees of freedom while E 2 depends on 
the second replica's degrees of freedom. In the classical limit, P{q) becomes the probability 
distribution for finding a fraction q of the degrees of freedom (x l ,pi) of the two replicas in the 
same elementary phase space cell. 

The operator inside the delta-function has a more concrete interpretation. For a given i, 

let 

rs rs 

Notice the r -B- s flip on the right hand side, which we can equivalently think of 

flip. Thus, what this operator X\ 2 does for a given i (when acting on a wave function say) is 

exchange the replica variables x\, x\. The first moment of P{q) is then 

i i i 

where S l (2) is the second Renyi entropy [91] of p l . If p = |\E')(\1/|, this gives information on 
how entangled the Xi are in the wave function if \1/ factorizes (no entanglement), then 
(q) = 1. For the isotropic ferromagnet (strong entanglement), we have (q) = |. The second 
moment is 

<« 2 > = w E i*<*xs>*i «-) = ]^E ^ w + ^ = ^ E e " s " ,2) • < 2 - 79 ) 

where p iJ is the density matrix obtained by tracing out all degrees of freedom except Xj, Xj, and 
in the last expression we dropped the 1/N term since it is understood that the thermodynamic 
limit iV — y oo is taken. The ratio (q 2 )/(q) 2 gives information about the degree of long range 
correlation (and failure of cluster decomposition) in the system. For a thermodynamic pure 



37 



state it equals 1, for the isotropic ferromagnet it equals 1 + Moments of order k are 
obtained analogously by fixing k degrees of freedom. 

Thus, for rather general systems, we have arrived at an infinite set of in principle com- 
putable order parameters that detect the emergence of a nontrivial state space landscape. 
Applications are discussed in [4, 5]. 

3 Supersymmetric quantum mechanics 
3.1 Introduction 

Complex systems are complex. As a result, it is hard in general to quantitatively determine 
even their basic features. When it can be done, such as for the Sherrington-Kirkpatrick model, 
the results often have a striking richness. Complex systems arising in string theory, such as the 
geometrically highly complex wrapped D-branes reviewed in section 4, are often significantly 
more intricate than simple spin models. This makes straightforward analysis more daunting. 
Although in principle the probability overlaps defined in section 2.5 define an intrinsic order 
parameter that could be used to analyze those stringy systems along the lines of the SK model, 
in practice this may still be difficult. This is even true for much simpler quantities like the 
entropy. 

On the other hand, the existence of dualities in string theory and the presence of supersym- 
metry often tremendously simplify things. In this section we will focus on supersymmetry. 
The main simplification that arises due to this is that many quantities of interest can be 
computed exactly in a semiclassical limit. The primary example of this is the Witten index 
[92], which counts ground state degeneracies in supersymmetric quantum mechanics, and is 
generically invariant under continuous deformations of the theory. Famously this allowed Stro- 
minger and Vafa [93] to microscopically compute the entropy of certain classes of extremal 
black holes arising in string theory, by tuning the string coupling constant from the strongly 
self-gravitating black hole regime to a weakly coupled, non-self-gravitating regime. In this 
regime the system has a well-understood D-brane description, and the Witten index can be 
computed in a relatively straightforward and precise way. 

The ground state degeneracy is not the only thing supersymmetry allows us to control. 
Through its relation with Morse theory, the structure of supersymmetric quantum mechanics 
gives useful insights into the saddle point structure of energy landscapes and the way they 
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are dynamically connected, and makes it possible to systematically compute nonperturbative 
level splittings and their associated exponentially slow relaxation rates. Clearly this is of 
potential interest in the study of complex systems. In fact, through a formal map between 
nonsupersymmetric classical statistical mechanics and supersymmetric quantum mechanics, 
this has ramifications beyond the realm of supersymmetric systems (see section 3.8). 

We will review these fundamental aspects of supersymmetric quantum mechanics in what 
follows. I will not assume any familiarity with the topic, and start from scratch, building up 
the topic along the lines of chapter 10 of [94], an excellent introduction to supersymmetric 
quantum mechanics and its interplay with geometry. A review focusing more on relations to 
exactly solvable systems is [95], and finally Witten's original [92] is a classic. 

Specific applications to string theory are left to section 4. 

3.2 Definition 

By definition, the Hamiltonian of a supersymmetric quantum mechanical system can be writ- 
ten as 

H = ±{Q,Q*}, (3.1) 
where Q, the supercharge, is an operator satisfying 

Q 2 =(g t ) 2 = 0. (3.2) 

The Hilbert space is Z 2 graded, splitting up in a even or "bosonic" part, and an odd or 
"fermionic" part: 22 

U = U B ®U F . (3.3) 

The operator (— 1) F is defined as +1 on %b and —1 on Hp- The operator Q is odd/fermionic, 
i.e. it maps from % B into %p an d from % F into 'Kb'- 

[(-1) F ,Q] = -Q. (3.4) 

Some immediate consequences are 

• [(— 1) F ,H] = 0, [Q,H] = 0, i.e. (— 1) F and Q are symmetries. 

22 At this level "bosonic" and "fermionic" are just conventional names. They do not necessarily mean that 
the quantum mechanical states represent bosonic or fermionic particles in the usual sense of the word, although 
in specific situations the notion may be correlated. 
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• H > 0, all energies are positive or zero. 

• H\a) = <^ Q\ot) = = Q^\a), zero energy states are supersymmetric. 

• For each positive energy level He — He,b © He,f, E > 0, the fermionic and bosonic 
subspaces are isomorphic. The isomorphism is given by the Hermitian supercharge 

Q 1 = Q + gt > (3.5) 

which satisfies Q\ = 2E on He and is therefore invertible when E ^ 0. 

The latter property is particularly important. It means that all positive energy eigenstates 
come in boson-fermion pairs. This is not necessarily true for zero energy ground states. The 
Witten index Q quantifies the degree to which this is not true: 

Q = dimHo,s — dimH ,F • (3.6) 

Since all positive energy states come in even-odd pairs, this can equally well be written in the 
following alternative ways: 

Q = Tr(-1) F (3.7) 
= Ti(-l) F e^ H (3.8) 
= J V(---)e- s{ - ) . (3.9) 

The last line represents the usual Euclidean path integral representation of the partition 
function, except that now both fermionic and bosonic variables satisfy periodic boundary 
conditions, instead of the usual antiperiodic boundary conditions on the fermions. This is due 
to the insertion of (— 1) F , which flips the sign of the fermionic variables. (The example below 
may make this more clear.) 

The most important property of the Witten index is that it is invariant under generic 
continuous deformations. This is clear from its definition: although the total number of 
supersymmetric ground states may vary when we vary the parameters of the model, as states 
move in and out of the E = level, all such arrivals and departures must come in bose-fermi 
pairs, and therefore the index will not be affected. 23 



40 



Figure 7: Id susy QM with h(x) = x 4 -5x 2 + 7. Blue: h(x), red: V(x) = ti(x) 2 /2, green: ® (x) cx e"^. 



3.3 One dimensional example 

A simple example is obtained by considering a bosonic (i.e. commuting) variable x and a 
fermionic (i.e. anticommuting) variable ■0, wave functions 

$(x,^) = &b{x) + $ f {x)$, (3.10) 

and canonical commutation relations [x,p] = i, = 1, where ip = ip*. So acting on wave 

functions we have 

p = -id X} ip = d^. (3.11) 

Then we consider the supercharge 

Q = tp(ip + h'(x)) , Q ] =ip(-ip + h'(x)) . (3.12) 

Here /i(x) is an arbitrary function, sometimes called the superpotential. Expressed in compo- 
nents $ = (* s ) this becomes 

Q=( ° °), (3.13) 



23 Although generically robust, there are important exceptions to this argument. In particular, when the 

gap we have implicitly assumed between the zero and first excited energy levels vanishes in the course of the 
deformation, the Witten index may jump. This can happen for example when the deformation passes through 

a degeneration allowing a ground state wave function to spread out all the way to infinite distance, becoming 

non-normalizable and therefore no longer a proper state in the Hilbert space, but rather part of the continuum 

of scattering states. 
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and the Hamiltonian 

H = \{Q,Q ] } 

= \{-d 2 x + h\xf + h"{x)[^M) (3-14) 

1 (-81 + h'(x) 2 - h"(x) 

2 y d 2 + h\x) 2 + h"(x) 

We can think of this as the Hamiltonian of a spin 1/2 particle of unit mass in a potential 
V(x) = \h\x) 2 and a magnetic field \h"(x). In the case h(x) = |x 2 , this becomes the 
harmonic oscillator potential V(x) = ^u 2 x 2 plus a constant magnetic field u. Assuming u > 0, 
the bosonic (i.e. $i? = 0) energy levels are then Eb = 0, u, 2oj, 3u>, . . . and the fermionic energy 
levels Ep = u,2u,3u, . . .. We explicitly see the pairing of E > states and the mismatch 
at E — 0, captured by the Witten index Q = +1. When on the other hand u < 0, we get 
E B = \u\,2\u\, . . . and E F — 0, \cj\, 2\cj\, . . ., so Q = — 1. 

More generally, it is easy to explicitly construct solutions to = Qt$ = 0: 

$ = Ab e~ h(x) + Ap e +h(x) $ 24 (3.15) 

For these to be true supersymmetric ground states, <3> has to be normalizable. When lim^i^oo h(x) = 
+oo, this requires A F = 0, and thus Q = +1. When lim^i^oo h(x) = — oo, we get instead 
Ab = and Q = —1. In all other cases, none of the solutions are normalizable, and Q = 0. 
An example is shown in fig. 7. 




3.4 More degrees of freedom 

We can easily generalize this to multiple degrees of freedom x 1 and ip 1 , I = 1, . . . , N, leading 
to 2 iV -component wave functions 

9@,x)= Yl OT 1 " " " (^ N ) bN , (3.16) 

6i,-,6 N =0,l 

and a fermion number quantum number F counting the number of ip 1 . We take the super- 
charge to be 

Q = J2^(Wi + dih(x)), (3.17) 

24 Hcre V> is, as in (3.10), a Grassman variable. Equivalently we could write $ = Ab er h ^ |0) + 
Ap e +htyX ^ V>|0), with |0) the "fermionic vacuum state", i.e. ?/>|0) = 0. 
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which satisfies Q 2 = 0. The Hamiltonian is then 




(3.18) 



It is now no longer possible in general to find explicit solutions to Q$ = Q^<& = 0, but when 
h(x) is a nondegenerate quadratic function, say 



7=1 

then the problem just factorizes into N copies of the Id harmonic oscillator situation we solved 
above, and thus the unique normalizable supersymmetric ground state is 



This has fermion number F equal to the number of negative eigenvalues of the Hessian of h 
at the critical point x — 0. This number \i is called the Morse index of the critical point. The 
Witten index is f2 = (— 1) M . 

Using this, we can compute the Witten index for any h(x) that has only isolated, non- 
degenerate critical points. 25 Indeed, since the Witten index is invariant under continuous 
deformations that do not change the asymptotics of h(x), we can consider the deformation 
h(x) — > Xh(x) and compute Q in the limit A — > oo. This limit corresponds to a semiclassical 
limit. The classical supersymmetric ground states are given by the critical points of h(x). 
The corresponding perturbative 26 semiclassical wave functions are to lowest order Gaussians 
centered at the critical points, of the form (3.20). All of these remain supersymmetric to 
all orders in perturbation theory. (Roughly this is because perturbation theory around one 
critical point does not know about the presence of other critical points, and if there is only 
one critical point, the Witten index is 1.) However, nonperturbative tunneling effects may lift 
this degeneracy, just like tunneling splits up the doubly degenerate perturbative ground state 
energy level of a double well potential in ordinary quantum mechanics. But because all lifting 

25 Such a function is called a Morse function, and it corresponds to the generic case. We will always assume 

h is a Morse functions in this section. 

26 More precisely, the perturbation theory for each critical point of h is a series in inverse powers of A 

and is obtained by changing variables from x to y defined by x 1 = x\ + y 1 j\f\ so H = \(Hq + \~ x l 2 Hi + • • ■ ) 

where Hq is the harmonic oscillator part of the Hamiltonian, obtained by truncating h at quadratic order in 

y, and Hi the part corresponding to the higher order terms. 




(3.19) 




(3.20) 



r.LU!<Q 
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Fi gure 8: Critical points of the vertical height function on a pretzel, with fermion number (Morse index) 
indicated. The Witten index is £1 = 3 - 9 + 2 = -4. 

must occur in bose-fermi pairs, it will not affect the total Witten index, so we can compute 
the index simply by adding the contributions of all critical points x± of h: 

tt= (-l)^- (3-21) 

a;*: h'(as*)=0 

Checking this for the Id example shown in fig. 7, we see that indeed there h has two minima 
and one maximum, so = +1 + 1 — 1 = +1, in agreement with our earlier result. 

When h is a harmonic function, as is the case when it is for example equal to the real 
part of a holomorphic function W(z) of variables z a = x 2a ~ l + ix 2a , a = 1, • • • ,n = N/2, all 
critical points have Morse index n, and the Witten index is equal to (— l) ra times the absolute 
number of critical points of W. This situation occurs naturally in the description of wrapped 
D-branes in Calabi-Yau compactifications of type II string theory. 

3.5 Curved spaces 

On a compact iV-dimensional Riemannian manifold M, we can naturally define a supersym- 
metric quantum mechanics with Hilbert space given by the space of complex valued differ- 
ential forms $ = $o + ^idx 1 + Qijdx 1 A dx J + • • • on M. The elementary 1-forms dx 1 get 
identified with the Grassmann coordinates $ . The inner product of two differential forms 
$(<*) = J2f =0 ®h-i F dx h A ■ ■ ■ A dx lF , a = 1, 2 is 

($W|$W>= / $(D A *$< 2 ) = ]T / d N x^g I ^...g I ^$M, lF ^l JF , (3.22) 
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where gu is the metric on M. The fermion number F is the form degree and the supercharge 
in the absence of a potential is the exterior derivative 

Q = d = dx*Vi , (3.23) 

where V/ is the covariant derivative. The Hamiltonian is the Laplace-Beltrami operator 

A = - (dS + Sd) . (3.24) 



^ = ^ 7 A, ^ = g IJ ^ JJ = g IJ LOl , (3.25) 



Comparing with our previous notation, we have, at the level of operators, the identifications 

d_ 

dip- 

implying in particular the anticommutation relations {ijj 1 , ip" '} = g IJ . More explicitly, using 
the relation between curvature and covariant derivative commutators, the Hamiltonian can 
be written as 

H = ~9 IJ ViVj + \Rijkl^¥^ K ^ L ■ (3-26) 

The supersymmetric ground states are the harmonic differential forms, and the Witten index 
is the alternating sum of the number of harmonic forms at each form degree, which is nothing 
but the Euler characteristic of M: 

N 

n = Y,(-iyP(M)= X (M). (3.27) 

F=0 

The Betti number If is also equal to the number of homologically independent (N—p)— cycles. 
On a compact manifold we have moreover b p = b p . An immediate consequence of this is 
that Q = when N is odd. For a Riemann surface of genus g, we have b° = b 2 = 1 and 
b 1 = 2g, so X = 2 - 2g. 

Adding a potential can be implemented by putting 

Q h = e~ h de h = d + dh. (3.28) 

This gives a potential term V = g IJ dihdjh in the Hamiltonian. Switching on a potential will 

not change the Witten index, since the spectrum is guaranteed to be discrete on a compact 

space. In fact, it will not even change the absolute number of supersymmetric ground states 

in each fermion number. To prove this one uses the very useful general fact 27 that the space 

27 This is not hard to show. First, for positive energy states, we have that if Q\a) — 0, then \a) = 
j^iQ^Q + QQ^)\ot) = Q (fsla)): so au positive energy states are trivial in cohomology. For zero energy 
states on the other hand we have Q\a) = = Q'\a), so if \a) — Q\f3), then (a\a) = (f3\Qi\a) = 0, i.e. |a) = 0. 
So each zero energy state is nontrivial in cohomology. 
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of supersymmetric ground states is isomorphic to the Q-cohomology: 

kerQ 

rLsusy — ~ • (o.ZJ) 

This is combined with the observation that there is an isomorphism between the cohomologies 
of Q = d and Qh = e~ h de h provided by the map \a) — > e~ h \a). 

An example of a Riemann surface, embedded in M. 3 is shown in fig. 8. This Riemann 
surface is known as a pretzel and acquires a natural metric from its embedding in IR 3 . For the 
function h we take the height along the vertical axis. There are 3 minima, 9 saddle points 
and 2 maxima, each of which gives rise in the large pretzel limit to a perturbative susy ground 
state, here of fermion number 0, 1 and 2 respectively. We know we can compute the Witten 
index in perturbation theory, yielding Q = 3 — 9 + 2 = —4, so the Riemann surface is of 
genus 3, which is confirmed by visual inspection. We know also that &o = 1 = 62, and thus 
bi = 6. Therefore, since the Betti numbers count the number of supersymmetric ground 
states at given fermion number, several of the perturbative ground states must be lifted by 
tunneling effects: only one linear combination of the F = and the F = 2 states can remain 
supersymmetric, and six linear combinations of the F = 1. The remaining three F = 1 states 
pair up with two F = and one F = 2 state and get lifted. The lifting will be exponentially 
small in the pretzel size (which sets the scale of h). 



3.6 Lagrangian formalism 

So far we have described everything in a Hamiltonian framework. Semiclassical tunneling 
amplitudes responsible for lifting perturbative ground states are however most easily computed 
in a Lagrangian framework. The Lagrangian framework and its associated path integral are 
also most convenient to efficiently obtain a differential geometric expression for the Witten 
index, as the integral over M of a differential form made out of the curvature, namely the 
Euler density. 

The Lagrangian is obtained from the Hamiltonian in the usual way. We are primarily 
interested in the Euclidean Lagrangian. Starting from (3.26) with the addition of the potential 
term, this is 

L = ^gux 1 ^ + ^g IJ d I hd J h 

+g IJ ^ I D T ^ J + l -R IJKL ^^ J ^ K ^ L + Vjdjh $ V • (3.30) 
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where r is the Euclidean time conjugate to H and x = D^ 1 = dtip 1 + Tj K d t x J ip K , 
V/Vj = <9/Vj — r^jV^-. The supersymmetries act on the fields as £(•••) = [eQ + eQ^, • • • ]: 

fa' = e^ 1 -eip 1 (3.31) 
<ty, 7 = e (-i / - Y l JK $ J $ K + g IJ djh) (3.32) 
^ = e^-r^V + ^j/i), (3.33) 

leaving the action invariant. 

As an application we derive a differential geometric expression for the Witten index, which 
is essentially the generalized Gauss-Bonnet formula [96]. The path integral representation of 
the Witten index is 

n = Tr( y -l) F e- fSH = j ' VxVi>V<i[)e-ti ' Ldt , (3.34) 

with periodic boundary conditions on all fields: <j>{0) = 0(0) for (f) = x,ip,ip. In the limit 
(3 — > 0, only the constant trajectories contribute to the path integral. Thus, taking h = and 
working in an orthonormal frame, the path integral reduces to the finite dimensional integral 

q = 1 f d N x d N ^ d N i) e -| W^-W* ; (3.35) 



where the normalization can be fixed for example by comparing path integral and canonical 
expressions for the free propagator in flat space (x\e~P R °\x). Using the symmetries of the 
Riemann tensor and the Grassmann integral representation of the Pfaffian of a matrix M, 
Pf M = fdip 1 --- dip N e -\^ M ^ 3 (with (Pf Mf = det M), this becomes 

Q= [ PiTZ, n I j = —R I JKL dx K Adx L . (3.36) 
Jm 47r 

This expression, which gives the Euler characteristic of M in terms of the integrated Euler 

density, is known as the generalized Gauss-Bonnet theorem. 

3.7 Inst ant ons and lifting 

By completing the squares, the bosonic part of the action (3.30) can be brought into the form 

fT/2 i 

S B = s(h(x f ) - h(xi)) + / dr - 9lJ (x 1 - sg IK d K h) (x J - sg JL d L h) , (3.37) 

J-T/2 1 

where s = ±1. The first term comes from integrating a total derivative, and x(—T/2) = Xi 
and x(T/2) = Xf are the start and end points of the trajectory. This implies the lower bound 
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Figure 9: Steepest ascent gradient flow lines for the function h(x,y) — A sin(27ra;) sin(27ry) on the 2-torus 
defined by the identifications a; ~ x + 1, y ~ y + 1. The critical points A, B are minima, G, H are maxima, 
C, D, E, F are saddle points. Big green arrows are paths which give a positive contribution to Q, big purple 
arrows give negative contributions. Little red arrows indicate wave function orientations. 

Sb > s(h(xf) — h(xi)) for s — ±1, hence Sb > \h(xf) — h(xi)\, which is saturated for steepest 
descent or ascent trajectories: 

x 1 = sg IJ djh(x) , s = sign(/t(x/) — h(xj)) . (3.38) 

These gradient flow trajectories (with the fermions put to zero) are automatically solutions 
to the Euclidean equations of motion, and moreover from the supersymmetry variations 
(3.32)-(3.33) it can be seen that ascending (s = +1) trajectories are annihilated by the Q- 
supersymmetry, while descending trajectories are annihilated by the Q^-supersymmetry. Such 
minimal action trajectories, in the limit T — > oo, are called instantons. 

Instantons play a key role in the computation of nonperturbative lifting of ground states. 
For concreteness consider the example shown in fig. 9, a two dimensional flat unit torus with 
superpotential h(x,y) = A sin(27ra;) sin(27ry). There are 8 critical points: 2 minima, 4 saddle 
points and 2 maxima, giving Vt = 2 — 4 + 2 = = x{T 2 )- The Betti numbers for T 2 are bo = 1, 
b\ = 2, 6 2 = 1, so again some of the perturbative ground states must be lifted by tunneling 
effects. The relevant tunneling trajectories responsible for this are precisely the instantons 
introduced above. 

The general story goes as follows. Denoting the perturbative ground states by lifting 
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means that the matrix elements of the supercharge reduced to the perturbative zero energy 
level, 

Qii = {i\Q\j), (3-39) 

although identically zero in perturbation theory, is not identically zero when nonperturbative 
effects are taken into account. Notice that since Q has fermion number 1, can only be 
nonzero if the fermion numbers of \i) and \j) statisfy the selection rule Fj = Fj + 1. The 
matrix elements can be given a path integral representation. It can be shown that thanks 
to supersymmetry, this path integral (generically) only gets contributions from instanton 
trajectories annihilated by Q, i.e. paths of steepest ascent, interpolating between critical 
points j — > i with F{ = Fj + 1. Moreover the Gaussian approximation to this contribution is 
exact, resulting in the very simple expression 

Qij = Yl n -y e ~ hl3 ' ( 3 - 40 ) 

where the sum is over steepest ascent flows from j to % satisfying Fj — Fj — 1. The prefactor 
n 7 = ±1 and hij is the instanton action, hy = \h(i) — h(J)\. This simplification in the path 
integral is referred to as localization. We will not derive this here; see [94] chapter 10 for a 
detailed pedagogical account. 

The sign of n 7 is determined by orientation considerations. This is best explained by 
going back to the torus example. First we need to choose an orientation for the perturbative 
ground states |z). In the figure the critical points are denoted by the letters A to G, and 
we will denote the perturbative ground states by these letters too. Then \A) and \B) are 
0-forms, |C), \D), \E), \F) are 1-forms, and \G), \H) are 2-forms. For example, to lowest order 
in perturbation theory, 

\A) = «e^ (2T)2( " 2+ " 2) , a = v / 2^A, u = x-\, v = y-\, (3.41) 
\C) = ae -!N 2 (« 2 +« 2 )4, u = x + y-l v = x -y + ±, (3.42) 
\G) = ae-^ 2n)2{u2+v2) duAdv, u = x - |, v = y-\. (3.43) 

Note that whereas the top and bottom forms can be given a natural orientation inherited from 
the orientation of the torus itself, this is not so for the 1-forms. We choose them such that at 
the critical point, the form is positive in the positive x-direction, i.e. the local coordinate u 
that is such that h ~ — u 2 + v 2 near the critical point is taken to increase when going to the 
right. In the figure, the orientation of the local u-coordinate is indicated by the red arrows. 
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The arrow is a local frame for the unstable direction. More generally, the choice of ordered 
standard local coordinates near each critical point determines a local frame for the unstable 
directions of h at each critical point. On the other hand when moving along an instanton 
trajectory, we can transport this frame from the initial point to the final point, and then add 
to it a new vector pointing in the direction of the path, to obtain a new frame for the unstable 
directions around the final point (which indeed generically has one more unstable direction 
than the initial point, in the direction of the arriving path). If this new frame has the same 
orientation as the frame we already had there, then n 7 = +1, and otherwise n 7 = —1. 

Using these rules, and noticing that in this example Ah = A for successive critical points, 
we get for example Q\A) = e~ x {\D) + \F) - \C) - \E)), Q\C) = e- x (\H) - \G)), Q\G) = 0. 
In matrix form, we thus compute Q, and H = \{Q, Q^} (writing zeros as blanks): 



/ 



Q 



-1 1 

1 -1 



-1 1 
1 -1 



\ 



-1 1 

-1 1 

1 -1 

1 -1 



V 



/l 

-1 



H = 2e 



-2A 



/ 



\ 



1 

-1 



V 



/ 



The space of true supersymmetric ground states is isomorphic to the cohomology ker Q/im Q of 
the matrix Q, which is relatively easily computed. This is called the Morse- Witten cohomology. 
Alternatively, we can find the actual supersymmetric ground states (within the approximation 
of working within the finite dimensional space of perturbative supersymmetric ground states, 
i.e. to leading order in e _A ), as the joint null space of Q and Q\ or equivalently the null space 
of H. The reduced Hamiltonian H furthermore gives the dynamical transition rates between 
different perturbative ground states, and diagonalizing it gives the energy level splittings. 

The conclusion in our example is that the supersymmetric ground states in this approxi- 
mation are given by the symmetric combinations \A) + \B), \C) + \F), \D) + \E), \G) + \H), 
consistent with the Betti numbers (1,2, 1). The four antisymmetric combinations orthogonal 
to those all have lifted energies equal to 4e~ 2A . 

Of course in more general examples, the splitting pattern and linear combinations will be 
less simple. But in principle at least there is a general algorithm to compute it in terms of 
simple, generically isolated, "atomic" (AF = 1) gradient flows, which is quite remarkable; in 
nonsupersymmetric quantum mechanics in dimensions greater than 1, doing something similar 
is typically very hard if not impossible. 
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For complex energy landscapes h(x) in high dimensional spaces, it will typically be more 
appropriate to compute statistical, thermodynamic quantities, but again the structure of 
supersymmetry will lead to significant simplifications. 

As a check on the above results, and also to get further insight into the nature of the 
approximation, observe that it is actually easy to write down the exact bottom and top 
fermion number ground state wave functions; they are: 



The approximation we made replaces this to lowest order in perturbation theory by an even 
superposition of Gaussians centered at the critical points. This is indeed an excellent approx- 
imation at large A. 

Note that in this limit perturbative ground states are exponentially close to being real 
ground states, and physically the number of perturbative ground states would therefore be a 
better measure for the ground state degeneracy than the Witten index. However, at generic, 
order 1 couplings, this is no longer the case, and one may reasonably expect all states that are 
not protected to be lifted to finite energies, making the Witten index an accurate, physically 
meaningful measure for the number of true ground states. On the other hand, as many 
examples in string theory have taught us, one must of course be careful that in the apparent 
strong coupling regime there isn't a dual weak coupling description hiding, which would again 
potentially produce many new near-exact ground states, not captured by the Witten index. 

This is the case in particular in thermodynamic, "large N" limits. For example the SK 
model and other glassy systems develop many pseudo ground states in the thermodynamic 
limit. Although not supersymmetric, there is an interesting map from classical statistical 
mechanics systems to supersymmetric quantum mechanics systems, showing this is not just 
an analogy, but really the same thing. To this we turn next. 

3.8 Supersymmetry, stochastic systems and glasses 

Supersymmetric quantum mechanics can be interpreted as a BRST description of a classical 
stochastic system [70, 97, 98, 71, 72]. This goes as follows. Consider a stochastic dynamical 
Langevin equation for a particle with position x(t) in a potential h(x): 



F = 0) = e 



-h 



F = dim) = e +h Yli*P I ■ 



(3.44) 




(3.45) 
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where r\ is a Gaussian random noise function with variance a, P{r]) oc e~^/ 2a (the noise can 
be thought of as thermal fluctuations with temperature proportional to a). Starting from 
some initial point X{ at time t — 0, the probability density of finding the configuration at a 
point Xf at time t is given by the path integral 

/d d 
VxVr)e-ti^ v2 5(-^- + h'(x) - rj) | det(— + h"{x)) \ , (3.46) 

appropriately normalized and with boundary conditions x(0) = Xi, x(t) = Xf. The combina- 
tion of delta function and determinant ensures that the path integration over x simply picks 
out with unit weight the unique solution to the first order linear differential equation (3.45) for 
given rj. It is the functional analog of the finite dimensional formula J dx 5(A-x+B) | det = 1 
where A is a constant matrix and B is a constant vector. Ignoring for the time being the abso- 
lute value signs, we can write the determinant as a fermionic path integral, so that amusingly, 
after also integrating out rj, we obtain exactly the Euclidean propagator of supersymmetric 
quantum mechanics with superpotential h(x): 

P(x f ,t\ Xi ,0) = J VxV^V^) e~ & ^ {±+h ' {x))2+ ^ + t> h " {x ^ (3.47) 
= e (h(xi)-h(x f ))/a J DxD^D^e-So ±(x 2 +h'(x) 2 )+4>i>+4>h"(x)i> ^ ^ 
= e (h(xi)-h(x f ))/ a ( X/ | e -"W/a| x .} _ (3_ 49 ) 

The last line is the propagator between states in the zero fermion number sector, with elapsed 
Euclidean time t and Planck's constant H = a (before we were setting h = 1). Inserting a 
complete set of (F = 0) energy eigenstates |n), we have 

(x f \e- tHs ^ /a \xi) = ^(x f \n)(n\xi)e- tEn/a . (3.50) 

n 

In the long time limit t — > oo, only the ground state survives. As we have seen, when 
Z = J dx e~ 2h ^l a < oo, the Witten index is 1 and the unique supersymmetric ground state 
is (x|0) = -^e~ h{x - la . We conclude 

P(x f , oo\x t , 0) = V^/) , p = -, (3.51) 
Z a 

with Z = J dxe -1311 ^. This is indeed the expected Boltzmann equilibrium distribution with 
inverse temperature /3, which alternatively can be derived from the Fokker-Planck equation 
associated to the Langevin equation. In fact the Fokker-Planck equation can easily be seen 



52 



to be equivalent to the Schrodinger equation in the zero fermion number sector (identifying 
P(x) oc e h ^ip(x)), providing an alternative derivation of the above identifications. 

Note that when the Witten index is —1, i.e. e +h ^ is normalizable, there is no future 
equilibrium state, but there is a past equilibrium state. When it is zero, there is no future 
and no past equilibrium. If there are multiple perturbative supersymmetric ground states 
with fermion number zero (i.e. local minima of h(x)), then in the small a (i.e. large /3) limit, 
there will be highly metastable states. According to (3.50), the time scale for decay of these 
metastable states is set by the nonperturbatively lifted energy, which as we have seen in the 
previous section can be computed from simple instantons and the Morse- Witten complex. 

All this generalizes in a straightforward way to higher dimensional and curved spaces. 
If h(x) is the Hamiltonian of some classical statistical mechanical system, then in the ther- 
modynamic limit N — > oo there may be degenerations that change the number of exact 
supersymmetric ground states at a given fermion number and even the Witten index. This is 
because degeneracy lifting instanton effects may vanish altogether, thus turning perturbative 
ground states into exact ground states. When this happens at fermion number zero, ergod- 
icity is broken and we get multiple exact equilibrium states. Whether or not this happens 
will depend on the distribution of energy barriers, i.e. the distribution of supersymmetric 
ground states with low fermion number, and how they are connected by instantons, i.e. the 
Morse- Witten complex. This should capture, if present, the ultrametric structure of the state 
space. 

Some discussion of spin glasses in the language of supersymmetric quantum mechanics 
can be found in e.g. [71, 99, 72, 59]. Recent considerations in string theory involving the 
above connection between stochastic equations and supersymmetry include [100, 101]. A 
nice example of the relation between transition rate distributions and universal properties 
of relaxation dynamics can be found in [102], where the typical log(l + t w /t) behavior of 
return to equilibrium after a perturbation of duration t w of a glass was derived from general 
arguments based on the form of a transition matrix with exponentially suppressed off-diagonal 
matrix elements. Such a transition matrix and Hamiltonian are typical for systems with a 
wide distribution of relevant instanton actions. 
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4 D-brane landscapes 



4.1 Introduction 

We now move on to the description of complex systems that appear in the context of string 
theory, and review some of the tools that have been used to describe their basic features, such 
as ground state degeneracies and distributions. One of the interesting things string theory 
adds to conventional analysis is the power of holographic duality [103], which in favorable 
circumstances can reduce a quantitative understanding of the strong coupling physics of the 
model under consideration to simple computations in gravity. This works particularly well for 
thermodynamics, as black holes are the holographic duals to thermodynamic states, capturing 
the strong coupling thermodynamics of the model with often rather astonishing simplicity. 

More specifically, we will have a look at the complex configuration spaces that arise when 
wrapping D-branes around various cycles of compactification manifolds, with various world- 
volume fields turned on. Such brane configurations arise in type II constructions of string 
vacua, where they are among other things responsible for the particle physics content of the 
compactification. They also arise in the description of charged black hole microstates, extrap- 
olated to weak coupling. In the weak coupling limit (large volumes, small string coupling), 
the only difference between the two as far as their perturbative description is concerned is 
the fact that in the former case, branes are filling the observable, noncompact space, while 
in the latter, they are localized at a point. But in particular the geometric internal space 
description of the configuration is mathematically identical in the two cases. Of course the 
actual dynamics and other physics of the two systems will be very different, one describing 
the universe we may live in, the other an object we may look at, but the fact that they are 
related at the level of classical configuration spaces is one of the beautiful incarnations of the 
unifying nature of string theory. Thus, in particular, the largeness of the landscape and the 
largeness of the entropy of a black hole are intimately connected. 

Despite this map, the fact that branes fill all of space has one dramatic consequence, and 
that is that the charges (brane wrapping numbers etc) are tightly constrained by so-called 
tadpole cancellation conditions. What this means is simply that whenever we have charged 
objects (such as a branes) with compact transversal dimensions (as is the case for space-filling 
branes but not space-localized ones), the sum of all charges must necessarily be zero. This is 
a consequence of Gauss' law, or more colloquially, the fact that flux lines otherwise have no 
place to go to; when there's a source, there must be a sink. Combined with the requirement of 
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stability (typically realized by staying close to a supersymmetric configuration), this imposes 
strict bounds on the allowed charges. 28 On the other hand, for space localized branes, no such 
charge restrictions arise, and in particular we can take the charges to infinity and consider 
thermodynamic limits, which can be expected to have universal properties. This motivates us 
to consider localized brane systems. 

At finite coupling another difference between the two creeps in: the presence of gravity. 
This has a very different effect in the two settings. In the space-filling case, "vacuum" solutions 
with slightly different internal brane configurations can give rise to very different asymptotic 
geometries. As emphasized by Tom Banks at this school [8], for various reasons including the 
fact that our most successful descriptions of quantum gravity use fixed spacetime asymptotics 
in an essential way, this makes it unclear if we can even think of these different solutions as 
being part of the same theory, or, as he would put it, it makes it clear that these solutions 
are not part of the same theory. This does not mean that the presence of a space of brane 
configurations is of no physical relevance here of course; they do all exist after all, and it is 
certainly of interest to classify the set of all possibilities. Furthermore local fluctuations as 
well as larger fluctuations in the form of bubble nucleations can certainly probe — albeit in a 
limited sense — the configuration space. But it does mean that the situation is conceptually 
much more subtle than, say, a collection of manganese atoms in copper. In contrast, for 
space-localized branes, these issues do not arise. Gravitational backreaction just turns the 
weakly coupled brane system into a black hole, leaving the asymptotics unaffected. Even 
better, thanks to the effects of gravity, we get a host of thermodynamic parameters for free, 
courtesy of the Bekenstein-Hawking entropy formula. This gives another motivation to study 
brane systems associated to black holes. 

Besides their intrinsic interest as models of complex, glassy systems with intricate land- 
scapes and known holographic dual descriptions, the study of such localized brane systems 
can also teach us things about the space-filling systems. Given their close physical and math- 

28 These bounds are set by the background curvature or orientifold contributions to the charge, which in 
turn are set by topological invariants such as the Euler number of the compactification manifold. Although 
these numbers can get fairly high, for example there are known [2] elliptic CY 4-folds suitable for F-theory 
compactifications allowing up to O(10 5 ) transversal D3-branes, constructions get increasingly sparse at the 
higher end, leading many to suspect the set is finite. For reasons we do not understand any deeper than 
through case by case examination of various conspiracies, the string theory landscape is uncannily "boxed 
in", making it for example very hard, and plausibly impossible [38], to find infinite families of parametrically 
controlled solutions with positive cosmological constant. 
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ematical relation, it is clear that they contain information about each other. For example, 
their Witten indices, suitably defined, are the same. Indeed, Witten indices of supersymmet- 
ric quantum field theories are traditionally computed by reduction to a quantum mechanical 
model. In the present context one can go further, and use the computational power of gravity 
(manifested through for example the Bekenstein-Hawking formula or its refinements [104]) to 
extract information about ensembles of possible string vacua. This idea was made precise 
in [105]; I will not have much to say about it in these lectures, except that it serves as yet 
another motivation to study space-localized branes. 

The analysis of D-brane systems as mean-field glassy systems, or in other words connecting 
this section and section 2, will be addressed in [5]. 



4.2 Low energy physics of the weakly coupled wrapped D4 

In this section we derive in detail the quantum mechanics describing the low energy dynamics 
of a D4-brane wrapped on a 4-cycle, in the limit of small string coupling and large volume. 
Besides being an interesting and accessible example of a complex system in string theory with 
applications to black hole physics, mathematics and string model building, it also provides 
a clean example illustrating various useful Kaluza-Klein reduction techniques, in particular 
the relation between fluxes and superpotentials. (The idea to consider such an open string 
landscape was proposed in [106].) 

4.2.1 Actions 

Consider type IIA string theory on IR 4 x X, where X is some compact 6-dimensional manifold. 
The part of the IIA bulk action relevant for our purposes is, in units with the string scale 
4 = 2iry/a' = 1: 

•S'liA = SliA,NSNS + <SnA,RR (4-1) 

Si,a,nsns = 2tt [ d 1G x^g(^- 2 R-j 2 H MNP H MNP \ (4.2) 



>IIA,RR 



JR 4 xX V z z " / 



(4.3) 



where H = dB is the NSNS 3-form field strength, F^ +2 ^ = dC^ p+1 ^ are the RR field strengths 
and g s is the string coupling constant. The theory contains Dp-branes with p even. In 
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particular a single D4 wrapped around a 4-cycle E in the compact manifold X has, in the 
limit of large volume and weak string coupling, an effective worldvolume action 



Sd4 — Sd4,nsns + Se>4,rr (4.4) 

Sd4,nsns = -— f d 5 xV^h(l + lFjjF IJ ) (4.5) 
9s Vmxe V ^ / 

5 D4 ,rr = 2n [ FAC& + (]f A F + ^-e(R)) AC* 1 ', (4.6) 

where hjj is the induced brane metric, F = dA + B, A is the U(l) gauge connection living 
on the brane and e(R) is a quadratic polynomial in the tangent and normal bundle curvature 
forms. When X is a Calabi-Yau manifold, this is the Euler density e(R): 

e{R) = Pin, TZ I j = ^-R I JKL dx K Adx L , (4.7) 

47T 

which we encountered earlier in (3.36). 

The couplings to the RR fields given by Se>4,rr are topological in the sense that they are 
metric-independent. They specify the gauge flux and curvature induced D2- and DO-charges 
on the brane. 

We have ignored curvature corrections to Se>4,nsns5 but in the end we will fix this by 
matching the energy required by supersymmetry for a given charge. A slightly more precise 
version of this action and references to derivations in the literature can be found in section 2 
of [107]. 



4.2.2 Bulk Kaluza-Klein reduction 

The low energy dynamics of the wrapped D4 system is obtained by Kaluza-Klein reduction. 
This is done most easily, and leads to the richest structure, when there is residual supersym- 
metry. An M = 1 theory in four dimensions does not have pointlike objects that preserve 
part of the supersymmetry, so one has to consider at least M = 2, i.e. compactifications with 
eight preserved supercharges, with the wrapped D-brane breaking half of those. When no bulk 
magnetic fluxes are turned on, this requires the compactification manifold X to be Calabi- 
Yau, i.e. Ricci-flat and Kahler. Concretely this means that there exist complex coordinates 
y m , m = 1, 2, 3 and a holomorphic complex orthonormal frame 8 m = 8 m n dy n , such that the 
metric is of the form 
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with associated Kahler form J and covariantly constant holomorphic 3-form Vt given by 

j = - (o 1 a e 1 + e 2 a e 2 + e 3 a e 3 ) , n = ^A^ 2 A^ 3 . (4.9) 

The Calabi-Yau condition means that J and Q are globally well-defined. The Kahler form 
has one holomorphic and one anti- holomorphic index and is therefore called a (1, l)-form. 
Similarly Q is a (3, 0)-form. With the above definitions we have 

i [ nAH = 8 f — = 8V X , (4.10) 
Jx Jx 6 

where Vx is the volume of X. 

As a concrete example we can take X to be the quintic Calabi-Yau, described by a homo- 
geneous degree 5 polynomial equation in CP 4 . The four complex dimensional projective space 
CP is the set of all nonzero (x\, 22, £3, £4, £5) € C 5 , modulo the equivalence relation x ~ Xx, 
A G C*. The simplest choice of polynomial is of Fermat form: 

X : Q 5 (x) = x\ + x 5 2 + x 5 3 + x\ + x\ = 0. (4.11) 

Deformations of the metric preserving Calabi-Yauness are called (geometric) moduli. The four 
dimensional low energy effective field theory will contain these moduli as massless scalars. In 
general it is not possible to write down explicit expressions for the Ricci flat Calabi-Yau metric, 
but Yau's theorem states that for a given Kahler class (i.e. the cohomology class of the Kahler 
form) and complex structure (i.e. a choice of complex coordinates), there is a unique Ricci flat 
Kahler metric. Thus, there are two kinds of moduli, those specifying the cohomology class 
of the Kahler form and those specifying the complex structure. In the case at hand there is 
just one Kahler modulus, which parametrizes the overall scale of the metric, i.e. the volume 
of X — the second cohomology of X is one-dimensional. The different complex structures 
are simply parametrized by the choice of defining polynomial Q$(x), modulo linear coordinate 
transformations of the Xj. This space is 101-dimensional. 

More generally, the dimension of the second cohomology of X will be 6 2 (X)-dimensional, 
i.e. there are b 2 (X) independent harmonic 2-forms on X. Equivalently, because of Poincare 
duality, 29 there are b 2 (X) independent homology 4-cycles in X. Let Da, A = 1,...,B 2 = 

29 Poincare duality is a natural isomorphism between p-cycle homology classes and (n — p)-form cohomology 
classes in an n-dimensional space. For ap-cycle C locally given by n—p equations fi(x) = 0, the associated dual 
cohomology class can be represented by the closed delta-function (n-p)-form C — 6(fi) dfiA- ■ -A<5(/ n _ p ) df n _ p . 
Note that for two cycles C and C with dimension adding up to n, we can thus write J C A C' = J C C' = 
#(CnC"), where the last expression counts intersection points for generic representatives, with signs according 
to orientations. Usually we will drop the hat in the dual to avoid cluttering. 
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b 2 (X) be an integral basis of harmonic 2-forms. Then 

J = J A D A , (4.12) 

and the J A are the Kahler moduli. We choose the signs of the 2-forms Da such that J A > 0. 

The reduction of the RR potentials produces, among other fields, electrostatic potentials 
in four dimensions coupling to D2 and DO charges: 

C (3) = <f) A dtAD A + --- , C (1) = <f)° dt + ■ ■ ■ . (4.13) 
4.2.3 Brane Kaluza-Klein reduction 

In the weak coupling limit g s — > 0, the backreaction of the wrapped D4 on the bulk geometry 
is negligible, due to the lower power in l/g s appearing in the D4-action. We can therefore 
consider the D-brane to be a probe in a background geometry specified by arbitrary, fixed bulk 
moduli. For the wrapped brane to be supersymmetric, it must wrap a minimal volume 4-cycle, 
or somewhat stronger even, it must be holomorphic [108]. In the case of the quintic this is 
equivalent to saying it can be described by a degree N homogeneous polynomial equation 

£ : P N (x 1 ,x 2 ,x 3 ,x 4 ,x 5 ) = 0. (4.14) 

The degree N can be identified with the multiplicity of the D4-charge. For example for 
N — 1, the most general polynomial is P\{x) = a\X\ + • • • + a^x^. The complex coefficients 
c are deformation moduli of S. Since overall rescaling of the coefficients does not change 
the zero set, we have a four complex dimensional deformation moduli space, topologically 
CP 4 . An example of a degree N polynomial is {P\{x)) N , which corresponds to N coincident 
branes wrapping Pi(x) = 0, with gauge group enhanced to U(N). The most general degree N 
homogeneous polynomial on the other hand gives a smooth singly wrapped brane with gauge 
group U(l). In the case at hand it has d(N) = + — 1 independent deformation moduli 
z a , as can be verified by direct monomial counting, taking into account that when N > 5, 
adding a term qN-b{x)Q^{x) to Pn(x), with Q 5 as in (4.11), does not change the zero set. 
More generally one can deduce the number of deformations using standard algebraic geometry 
techniques [109]. For a D4-brane wrapping a homology class 30 [E] = N A D& with N A > this 
gives 

d = ]:D ABC N A N B N C + ^-c 2>A N A - 1 . (4.15) 
o 12 

30 Here the Da are the Poincare duals to the Da introduced in (4.12), for which as mentioned in footnote 
29 we use the same notation. 



59 



Here the "triple intersection numbers" Dabc are defined as 

D A bc= [ D A AD B AD c = #{D A nD B nD c ) (4.16) 



x 



and the C2,a are topological numbers (the second Chern class of X) depending on X only. 
For the quintic we see by comparing that D m = 5 and = 50. The same techniques allow 
straightforward computation of the generic Euler characteristic of E, namely 

X = D ABC N A N B N C + c 2A N A , (4.17) 

which for the quintic leads to % = 5 N 3 + 50 N. The second Betti number, i.e. the number of 
independent 2-cycles or 2-forms on E is obtained from this as 62(E) = x ~ 2 — b\ — 63. When 
N A > the Betti numbers b\ = 63 of E are inherited from X, so they vanish except in the 
higher susy cases when X is T 6 (b\ = 6) or T 2 x K3 (b\ =2). In any case, 62 grows as iV 3 , and 
therefore, crucially, we will get a huge magnetic flux degeneracy on these branes in the large 
N limit. Indeed we can turn on harmonic gauge field strengths F = dA on E characterized 
by 62 integers Si, the flux quanta: 

62(E) 

F=Y,S l( r\ (4.18) 

i=l 

where the <j % form a basis of integrally quantized harmonic 2-forms on E. In terms of Poincare 
dual 2-cycles, we can equivalently write 



Si = rfo If A o* = Vij [ F, V ^=[a l A ^ = #K n a J 
J Jul Jt, 



(4.19) 



with 7]ij the inverse of rfi . The matrix rfi is called the intersection form of E. It is integral 
and unimodular, but in general not positive definite — in fact as we will see later at large N 
it has signature (b^, b%) oc (|, |) x iV 3 . 

The flux and curvature induced D2 and DO brane charge can be read off from (4.6) with 
the reduction (4.13): 



5 D 4,RR = ^ J dt (-qo<p + q A <p A ) 



(4.20) 



— X(E)-- / F A F = —— — - 
24 AV ; 2 7 S 24 2 

q A = I FAD A = D\Si , (4.22) 



y. 
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Z{z + 5z) 



Figure 10: Infinitesimal variation of £ along its moduli space. 

where the integers D\ are given by 

D A ,i = J <r' A D A = #K n D A ) . (4.23) 

Kaluza-Klein reduction of the D4 produces a supersymmetric quantum mechanics with 
c£(E) complex continuous degrees of freedom z a and 62(E) discrete flux degrees of freedom 
The moduli z a parametrize a moduli space .M. The discrete fluxes can be thought 
of as quantized momenta dual to periodic coordinates that do not explicitly appear in the 
Hamiltonian. The wrapped D-brane with F = preserves 4 supercharges; it is the dimensional 
reduction from 4 to 1 of a four dimensional M = 1 theory In a sector with fixed flux we 
therefore expect the low energy effective action for the z a to be of the general chiral multiplet 
form [110] 31 

S=—Jdt\ [g a iz a f - g ai d a W(z)B- b W(z)) +■■■, (4.24) 

where the + • • • consists of terms independent of z, W(z) is a holomorphic superpotential 
(absent when F = but generically nonzero when F^O), and g a i a Kahler metric: 

g ai = d a B- b K(z,z). (4.25) 

We have explicitly retained the overall prefactor — instead of absorbing it in the metric and 
superpotential. In this way g s can be thought of as Planck's constant in the supersymmetric 
quantum mechanics. 

31 Note that this is of the form (3.30) (in its Lorentzian version) with h = \{W + W). The special form of 
h is due to the presence of four supersymmetries Q a , Q' a , a = 1,2, rather than the generic minimal amount 
of 2 supercharges. 
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To identify the metric, we should take the brane to be slowly moving along its moduli 
space and expand the first term in (4.5) to second order in the velocities. An infinitesimal 
displacement 5z a along the moduli space causes an infinitesimal normal displacement 5y m = 
v™8z a of E inside X, where the v™ depend holomorphically on the coordinates of E. This is 
illustrated in fig. 10. Thus, when slowly moving along the moduli space with velocity z a , the 
first term in (4.5) becomes 

- / d 5 xV^h= I dt I dv(-l + y m nV™vfz a A , (4.26) 

where dV is the volume element on S. Choosing the orthonormal frame 9 appearing in (4.9) 
such that 9 1 and 9 2 lie along E (i.e. they span T*E) while 9 3 is normal to it, we can write 
dV = (l^ 1 A 9 1 ) A [\9 2 A 9 2 ) so the kinetic term becomes 

gl A /)2„,3 a fll a fl2— 3 | z,a£b ( I , a _ , \ " 

S 



A 9 2 vi A 9 A 9 2 vl z a t b = - / w„ A w d i a l b , (4.27) 



where 

uj a = VL-v a = 9 1 A9 2 vl (4.28) 

is the contraction of Q with the vector field v a . It is a well-defined holomorphic (2, 0) form on 
S. This map is an isomorphism between deformations of S and harmonic (2, 0) -forms on E. 
By comparing to (4.24) we thus identify 

9al = 7 / w a A Q~ b . (4.29) 



4 „ 

To show that this is Kahler, we first expand the u a in the integral basis introduced in (4.18): 



U a = 0J ai (J % , UJ ai = Tjij / 0J a ■ (4.30) 



Now fix an arbitrary reference point z a = in the moduli space M. parametrizing the holo- 
morphic cycles S(z). If we move away from this point, the 2-cycles a % will sweep out 3-chains 
P(z). This is illustrated in figure 11. We can integrate the holomorphic 3-form over these 
3-chains, producing the holomorphic 3-chain periods 



W(z) = Q. (4.31) 

JY^z) 

Because E is holomorphic this is independent of the choice of representative of a 1 inside E(z). 
By the definition of the u a , we then have 



d a W{z) = / n-v a (z) = / u a , (4.32) 

JlJ i {z) J (ji 
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Figure 11: The inner (yellow) hyperbole-id represents S(0), the outer (blue) one is £(z), the green horizontal 
disk stretched between the two represents the 3-chain T(z), and its boundary in E(z) is a(z). 



and combining this with (4.30) and (4.29) gives us 



-daU'^diW =d a d- b K(z,z) 



1 



ITTfalF. 



(4.33) 



9ab — 



K 



4 



This makes it manifest that the metric is Kahler. 

To identify the superpotential, we consider the second term in (4.5), which for a given 
magnetic flux F gives a potential energy 



where * is the Hodge star operator on E and we used * 2 = 1. The second term is topological 
- it does not depend on the z a , and (4.21) shows it is nothing but the DO-charge induced 
by the fluxes. Thus, the magnetic potential energy is bounded below by the flux induced 
DO-charge, with the bound saturated by anti-self-dual flux F = — * F . In the following we 
will show that the first term is (almost) the |cW| 2 term of (4.24). 

The -k operator commutes with the decomposition of 2-forms according to their number of 
(holomorphic, antiholomorphic) indices. 32 It acts as +1 on (2, 0) forms: This 
is easily checked by considering a holomorphic orthonormal frame of X: -kiO 1 A 9 2 ) = +6 1 A 9 2 . 
The complex conjugate (0, 2)-forms are self-dual as well. For (1, l)-forms there are two cases to 

32 This is true because E is a Kahler manifold, being a complex submanifold of the Kahler manifold X. 




F A F , 



(4.34) 
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consider: (1, l)-forms proportional to the Kahler form are again self-dual, while all (1, 1) forms 
orthogonal to that are anti-selfdual, i.e. * = — 1. This is checked by considering 6 1 A6 1 + 8 2 A0 2 

resp. e 1 a e 1 - e 2 a e 2 . 

Accordingly the first term in (4.34) is 

f ](F + *F) A(F + *F) = I 2F™ A F^ + F^ 1] A f[ 1a) . (4.35) 
Jt, 4 Jt, 

The component F+'^ is obtained by projection onto the Kahler form J = J A D A (pulled back 
to E): 

Fi 1A) = f + J, f+ = J—^FAJ=jJ-^q A J A , ^J 2 = D ABC N A J B J C , (4.36) 

where we used (4.22). This expression is again independent of the D4 moduli z a . The second 
term in (4.35) thus equals -j^-p (q A J A ) 2 - Note this is independent of the overall scale of J. 

The remaining term does depend on the moduli z a and we will see it can be identified with 
\dW\ 2 . Expand F^ in the basis of (2, 0)-forms defined in (4.28): 

F^ = f a u ai f a = \g al J FAur b , (4.37) 

with g a i defined in (4.29) and g ah its inverse. Using the decomposition (4.18) and the relation 
to chain period derivatives (4.32), we write 

f FAur b = Si f u~ b = Bb {S t W(z)) , (4.38) 

and the first term in (4.35) becomes \g ab d a (SiTl % )di(SjTV). Comparing to (4.24), we identify 
the flux induced superpotential 

W(z) = SiU^z) . (4.39) 

We put everything together below. 

4.3 D4 supersymmetric quantum mechanics 
4.3.1 Bosonic action 

To summarize, the low energy degrees of freedom at weak coupling and large volume of a 
generic 33 D4-brane wrapped on a 4-cycle £ = N A D A in a generic Calabi-Yau 3-fold with 

33 Generic means here that the brane world volume is smooth and the gauge group is U{\). Coincident branes 
will give rise to an enhanced nonabelian gauge symmetry and intersecting branes will give rise to additional 
light bifundamentals. We also do not include pointlike DO-branes yet. We will return to this in section 4.4. 
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triple intersection numbers D AB c an d Kahler form J = J A D A (J A , N A > 0) are d = b 2+ — 1 
complex deformation moduli z a and b 2 = X ~ 2 integers Si parametrizing the U(l) flux 
F = S i( j\ with b 2+ = \D ABC N A N B N C + ±c 2A N A and X = D ABC N A N B N C + c 2)A N A . Up 
to a constant energy term the bosonic part of the action is 



S = S dyn + S top (4.40) 

s dyn = — \g*?? - \g al d a wB- b w (4.4i) 

^to P = -—Jdt[(q + (q A J A ) 2 ) + g s (go 0° - QA <P A )] , (4.42) 
with superpotential and Kahler metric given by 



W(z) = SiW(z), W(z) = / Q (4.43) 

JV i {z) 

g al = d a d- b {W VijW) , (4-44) 



flux-dependent DO- and D2-charges 



» = -i-h" s - s ' = -Yi-\l FAF (445) 



s 



q A = D* A St= I D A AF, (4.46) 

and unit Kahler form 



J (D ABC N^jBjCy/2 J (J s J 2) l/2 J - ( 447 ) 

The matrix rfi = ^(crTlff') is the intersection form of S and .D^ = #(D A r\a l ). The 3-chains 
T*(z) are 3-chains swept out by the 2-cycles a 1 on E when moving from a fixed reference point 
to z in the moduli space. We recall that for the quintic D 1U = 5, J 1 = c 2> i = 50. 

This is a reliable description at weak coupling and low energies. What exactly does weak 
coupling mean? We recall from our general discussion of supersymmetric quantum mechanics 
in section 3 that weak coupling corresponds to a limit in which the superpotential h is scaled 
up as h — > Xh with A — > 00, while keeping the metric constant. To see how this relates 
to other scalings, consider the simplest case of a single variable with bosonic Lagrangian 
L = Ax 2 + B(dh/dx) 2 . Then we can redefine x = y/y/A such that L = y 2 + AB(dh/dy) 2 . 
So weak coupling means AB — > 00. Applying this to the case at hand and remembering that 
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because of (4.10) the metric and periods scale as g a i oc Vx, IP oc y/Vx with the volume of X, 
we have effectively A oc Vx/g s and B oc l/g s - Thus weak coupling means 

Vx 

— TP — y oo . (4.48) 

9s 

The left hand side is essentially the ratio of the string length squared over the 4d Planck length 
squared, t? s /£\- Hence the weak coupling regime is the regime in which the string length is 
much larger than the 4d Planck length, and the low energy supergravity description breaks 
down. 

We see from (4.41) that switching on the flux quanta Si generates a highly complex po- 
tential energy function for the z a , leading at fixed go and g^ (and large go) to a vast energy 
landscape with exponentially many minima, parametrized by the values of Si and z a for which 
dWs(z) = 0. This is kinematically very similar to the landscape of flux vacua in string theory, 
except of course that here we have a quantum mechanical system rather than a universe. 

The potential is in general not a single valued function on the moduli space Ai, due to 
possible monodromies a 1 — > M^o - - 7 of the 2-cycle basis when going around noncontractible 
loops in Ai. It does become single valued on the covering space Ai, also known as Teichmiiller 
space, but then one has to quotient the theory by the covering group, which acts as a discrete 
gauge symmetry on z a and Si. The metric g a i and the charges go and g^ on the other hand 
are single valued, since they can be expressed without reference to the basis er\ 

The different zero energy classical minima of the potential can mix quantum mechanically 
through tunneling instantons. By completing the squares in the Euclidean version of the 
action one finds a bound S^ yn > y- Re(e~ ta AW) for any real a, saturated when 

z a = - e - ia g ai B- b W(z) . (4.49) 

The strongest bound is obtained with a = arg(AW / ), so this is what we should take a to be 
in order to be able to find instanton solutions. The instanton action is then 

2-7T 

Sg n = -\AW\, (4.50) 
and the trajectories are straight lines when projected to the H^-plane. 

4.3.2 Geometry of moduli space 

Before turning to the supersymmetric completion of the model, it may be useful to quickly 
review some basic formulae in Kahler geometry as well as more specialized expressions ap- 
plicable to the model under consideration. When a metric is Kahler, i.e. g a i = d a diK, all 
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Christoffel symbols Yfj = \g KL {QiQlj + &j9il — Ql9ij) vanish except when all indices are 
holomorphic or all indices are anti-holomorphic. When all indices are holomorphic, we have 

T a bc = g a %g cd - = |<9°ff = \f u a A d b co c . (4.51) 

Jx 

We are lowering and raising indices with g a i and 77^'. The first expression is true for general 
Kahler manifolds, the remainder specializes to (4.44). Equivalently, J ' u a A Vj,c<; c = 0. The 
Riemann curvature, defined by [V/, V j)X K = Rjjl K X l , also simplifies: 

Raid = QedhKc = \ v-^-rr v a <9 c ir = \ [ %u s a v a u c . (4.52) 

Jx 

Again, the first equation is valid for general Kahler manifolds. Components of the curvature 
tensor not related by symmetries to the above ones (such as R a bai) & U vanish. Finally, in 
addition to the usual curvature symmetries Rijkl = —Rjikl = —Rijlk = Rklij and 
Ri[jkl] = 0, a Kahler manifold also satisfies R a i cd - = R c b a d- 

Geometrically, the fact that J u; a A VfeW c = means that VftCo> c is of type (1, 1). Indeed, 
taking a derivative of a (2, 0)-form turns it into a (2, 0)-form plus a (l,l)-form (because the 
pullback of a (2, 0) form from E(z) to H(z+5z) will at most produce one extra antiholomorphic 
leg to first order in 5z). The same is then true for VfeO; c , and J o) a A VbOJ c = further implies 
that the (2, 0) part is actually zero, leaving only a (1, 1) part. 

Keeping in mind the (p, q)-type of various forms, the following orthogonality properties 
hold: 

,9JP<9 6 IL = 0, <9JT Vm9 c IT = , d a WD Ai = 0. (4.53) 

Further orthogonality properties can be derived from these by differentiation, for example 
V b d a W D Ai = 0. 

4.3.3 Hamiltonian and supersymmetric completion 

To find the supersymmetric completion of this model, we could either reduce the fermionic 
part of the full D4-action, or we can infer it from the structure of the bosonic part. We will 
go for the latter route. To do so, we first write down the bosonic Hamiltonian derived from 
the action given earlier. This is = H^ n + with 

<L = ^ a Wi + ^/m^5 3) (4.54) 

H£l = g(go + (J~W). (4-55) 
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Here p a = ^g^z 11 and we dropped the "chemical potential terms" proportional to 0° and 
<f> A here. If we consider the Si as constants, this is of the form of a 4d Af = 1 theory of 
chiral superfields reduced to Id, and the supersymmetric completion would be immediate. 
However the Si are not dynamically conserved in this system, since motion from a point z in 
the moduli space Ai to itself along a noncontractible loop will in general change Si — >■ M£ Sj 
due to monodromy. On the other hand go, Qa and H^ p are monodromy invariant, and hence 
they are conserved. In particular for an instanton, the change AW that appears in (4.50) is 
due to a change Az with constant Si on Teichmuller space, but it will in general lead to a 
change AS" when reduced back to the moduli space. 34 Therefore we should consider the Si to 
be dynamical. We can think of the Si as quantized momenta conjugate to angular coordinates 
d l that do not appear explicitly in the Hamiltonian, i.e. Si = —id$i. The d 1 can then be 
viewed as the potentials obtained by KK reduction of the 2-form potential on the D4 that is 
dual to the U{1) gauge field. The appropriate metric for these coordinates can be read off 
from the above expression for the Hamiltonian, which can be cast in the form 

= ^g a ~ b p a p h +^S l S,, (4.56) 

where we have defined a metric on flux space 

g lj = 2 (niTl ja + ntn^ + JV J ) - rf j , Ul = ^djV , ~f = J A D\ . (4.57) 

In this expression we use rjij and g a i for index raising and lowering. The three terms inside 
the brackets are projectors onto (2, 0)-forms, (0,2)-forms and the (1,1) direction parallel to 
the Kahler form; together they project to the space of self-dual forms. That these terms are 
properly normalized projectors follows from the relation (4.29) and from the definition (4.47). 

In fact this metric is nothing but the the Hodge ^-product, and the decomposition into 
projections is just (4.34) with (4.35) again. Consistent with this, we have that g l jg\ = 5 % k , 
which is the property * 2 = 1. This also shows that the inverse of is g^ = rf k rf l gki- 

The above Hamiltonian is just that of a free supersymmetric quantum mechanics in a 
curved space. Accordingly we could construct the usual supercharges Q = d and = 
d) and take H = \{Q-,Q^}-, i.e. the Laplacian. However then the Q = supersymmetric 

34 Incidentally, an open D2-instanton wrapping a special Lagrangian 3-chain T in X with boundary dT = 
<t + — cr_ on S has an action 5d2 = fr I Jr^l = f 5 " l^^l' wnere AW is the change in superpotential when 
F — > F + er + — er_ . This is exactly the same as for monodromy instantons, suggesting a possible identification 
of the two. 
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ground states of the system would have H dyn = H top = 0, which corresponds to D-brane 
states without any flux. This is too restrictive. We only need ifdyn = for a state to be 
supersymmetric, so we need to find supercharges that square to i^dyn only. These are basically 
the dimensionally reduced 4d Af — 1 supersymmetry generators, except that we interpret the 
Si now as dynamical momenta: Si = —idi. Introducing the fermionic operators ip a = dz a A, 
ip a = dz a A and their conjugates ip a = (ip a y , %p a = (ip a y , we define the supercharges 

Q- = ^r^a + y/fj 5 d- a fl l , Q + = ^rV- a + JfJ a dj? d t . (4.58) 

and their conjugates Q_ = (Q-)\ Q+ = (Q+V ■ Equivalently 

Q- = \[^ d + \[i mid n Q+ = \f^ B + ^ dul ^- ( 4 - 59 ) 

The supercharges satisfy the extended supersymmetry algebra 35 

{Qa,Qp} = 6 aP H dya , {Q a , Qp} = = {Q a , Qp) , (4.60) 

where i^dyn reduces to (4.54) in the zero fermion number sector. This defines the summetric 
completion of -Hdyn- Working out e.g. {Q_ , } we get explicitly, acting on wave functions 
<£>(?/>, a;) (i.e. differential forms): 

#dyn = -^fV a v a + ^R alcA r^ h ^ - f s d a wB a wd l d j - r^ b v a d b w d t + f b rv- a B- b Ti% . 

Here, again acting on differential forms, V a = d a + T b ac ip c ilJb, and p a = —id a , Si = —idi. 
Evidently neither the d l nor their fermionic superpartners (i.e. something like dd l ) appear 
explicitly in the Hamiltonian. Nevertheless they do not completely decouple from the z- 
degrees of freedom, again because of global monodromies — the ^-torus fibration over the 
moduli space M. has nontrivial identifications when going around nontrivial loops. 

To summarize, the full supersymmetric Hamiltonian including the RR potential terms is, 
in terms of the conserved charges: 

H = H dyn + 2 f s (g + ( J A q A f) + 2vr(0 o g o + <Aa) • (4.61) 

35 Because a 4d Af = 1 theory has four supersymmetries, the system has twice as much supersymmetry as 
the minimal susy QM case studied in section 3. As "the" Q-operator (satisfying Q 2 — 0, {Q, Q^} — H^ yn ), we 
can therefore choose from a family of linear combinations of the supercharges. The phase a appearing in the 
instanton flow equation (4.49) determines which supersymmetry is preserved by the instanton, and therefore 
which supersymmetry gets corrected. 
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Making the § momenta Si explicit in the topological part of the Hamiltonian, this becomes 

H = H dyn -2n(j- s + 0°) (* + l^SiSj) + f a (TStf + 27T(f> i S i . (4.62) 
Unpacking the whole thing gives us 



- <t>°V ij )SiSj + 2nL l S t . (4.63) 

Here <p l = <f> A D l A , g l i is the metric (4.57), and 

V = l - ^-rij b V a d b W + ^-f b ^V- a B- b W . (4.64) 

It is often convenient to split quantities like Si and Li into self-dual and antiself-dual parts, 
as this leads to more transparent "left-moving" and "right-moving" expressions. Thus, for 
instance 

(j- 3 9 ij ~ = (i - 4>°)Sl - (i + 0°)S! , (4.65) 

where squares denote contraction with rj, i.e. X 2 = XiX % = rf^XiXj (so X\ > 0, X 2 < 0). 
The projectors introduced in (4.57) can of course be used to give an explicit expression for 
I. and hence also for X_ = X — X>. 



4.3.4 Supersymmetric Lagrangian 



Legendre transforming the Hamiltonian, with z a = f^, d 1 = jNjr, we get the corresponding 



Lagrangian 

9s 7T 3s 24 

where Af = ^ a + z b T a bc ip c . Note that L + = <p + as {V a d b U) + = 0. 

4.4 Constructing supersymmetric flux configurations 

We have constructed the supersymmetric quantum mechanics describing the low energy physics 
of smooth D4-branes with £7(1) gauge symmetry and abelian fluxes F turned on. In the ther- 
modynamic limit N — > oo we get 0(N 3 ) discrete flux degrees of freedom Si interacting with 
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Figure 12: Analog in 2d of the problem of fitting a holomorphic 4-cycle to contain a given collection of rigid 
holomorphic 2-spheres in a 6-dimensional space: here we are fitting a 1-cycle given by a degree 5 polynomial 
in x and y to contain 36 randomly generated points, thus fixing all moduli of the 1-cycle. 

each other in a complicated, nonlocal way through coupling to 0(N 3 ) continuous degrees 
of freedom z a . The effective Hamiltonian (4.56) is somewhat reminiscent of models for spin 
glasses or the Hopfield model for neural networks, but more intricate, and a priori without 
quenched disorder, making it more like a structural glass than a spin glass. 36 In any case this 
suggests the presence of many effectively stable states. We will now make this quantitative, 
ignoring for the time being the possibility of brane degenerations and nonabelian gauge group 
enhancement. 

Classical supersymmetric brane configurations are those for which H^Si = 0, i.e. dWp{z) = 
0. As the 3-chain periods are very hard to compute explicitly, it might seem hopeless to try 
to construct any of those configurations explicitly. But a little thought shows otherwise. The 
condition dWp = is equivalent to 

In other words if we switch on a flux F = SiU 1 , the brane moduli must adapt themselves 
to make the flux of type (1,1). A geometric interpretation of this is that the flux becomes 

36 A spin glass has quenched disorder built into the Hamiltonian. A structural glass has a simple Hamiltonian, 
and spontaneously generates its own disorder. This includes the material normal people refer to when they 
use the word glass. 
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Poincare dual to a linear combination of holomorphic 2-cycles. Indeed the Poincare dual 
of a holomorphic 2-cycle a is clearly of type (1,1) since for every (2, 0)-form u we have 
j a A u = j^u = 0. The last equation is true because a is the image of some map y m (u) 
with u a single complex coordinate, and so the pullback of a (2, 0) form will be proportional 
to du A du = 0. Thus, any linear combination (with arbitrary signs) of fluxes Poincare dual 
to holomorphic 2-cycle will be of type (1, 1). The converse is also true as a special (proven) 
case of the Hodge conjecture. The subspace of moduli space where extra integral harmonic 
forms become of type (1, 1) is sometimes called the (generalized) Noether-Lefschetz locus in 
algebraic geometry. 

This interpretation gives a relatively simple way of explicitly constructing classical super- 
symmetric ground states. Let us illustrate this with an example. Consider again the quintic 
Calabi-Yau Q$(x) = and let us assume that we have chosen the polynomial Q§ such that it 
contains the following holomorphic 2-cycle, parametrized by (u, v ) e CP 1 : 

7 : (xx,x 2 , x 3 , x 4 , x 5 ) = (u, -u, v, -v, 0) . (4.68) 

The Fermat quintic Yli x i = f° r example contains 7. The topology of 7 is that of CP 1 , 
i.e. it is a 2-sphere. Such "degree 1 rational curves" are generically isolated in a Calabi-Yau. 
This can be checked in the case at hand by simple counting: the most general linear map 
(u,v) — > (x\, ■ ■ ■ 5X5), has 10 coefficients, 4 of those can be eliminated by linear coordinate 
transformations of (u,v), and requiring Qq(x(u, v)) = for all (u, v ) produces 6 equations on 
the remaining coefficients (namely the coefficients of u 5 , u A v, . . . , v 5 after expanding out), 
so we generically expect a discrete solution set. Indeed in the generic quintic there are 2875 
isolated degree 1 rational curves. 

Now consider our single wrapped D4 (N = 1), described by a linear equation 

S:^aiXi = 0. (4.69) 

i 

Requiring S to contain 7 produces the equations ai = a 2 , a 3 = a 4 . Therefore, if we switch 
on a flux Poincare dual to 7 on S, rigidity of the rational curve means that in order to 
preserve supersymmetry, these equations must remain satisfied when deforming S. Hence 
the D4 moduli space is reduced from CP 4 to CP 2 . If we also require say the rational curve 
7' : x = (0, u, — it, v, — v) to lie on E, enforced by switching on F = 7', we end up with 
a i — 02 = a 3 = (14 = 05, i.e. all moduli have been lifted. A simple 2d analog of this is 
illustrated in fig. 12. By combining different pairs of the 2875 degree 1 rational curves we can 
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thus explicitly construct * = 4131375 different isolated supersymmetric configurations 
(assuming there are no duplicates or degenerates). In fact for each choice of pair (7,7') there 
is an infinite number of different flux states, obtained by switching on a more general linear 
combination of the harmonic forms 7 and 7': 37 

p = k 7 + k'y, (4.70) 

where k, k' G Z. Of course these are not all degenerate in total energy. They all have i^dyn = 
but they will have different H top . To compute this we should find the D2- and DO-charges. 
They are 

qi (F) = J^D, A F = k + k' , q (F) = -Y A -\S/^ F = ~2A + \ {k2jr ^ ' ^ 

with x — 55. Here we used the fact that the self-intersection number of 7 inside £ equals —3. 
This can be computed either by explicitly deforming 7 inside £ (necessarily nonholomorphi- 
cally) and counting intersections, or by using some algebraic geometry (as in section 4.3 of 
[107]). The total energy H = H top is given by 

= q + — = -— + -(k 2 + k' 2 ) + -(k+ k'f . (4.72) 
2tt yo 5 24 2 v ; 5 V ; V ; 

Thus we get a nondegenerate lattice of vacua for each generic pair of degree 1 rational curves. 

This can be generalized in various ways. Instead of degree 1 rational curves, we can consider 
degree d rational curves. The same simple counting argument as before indicates again that 
they are generically isolated: there are 5(d+l) — 4 coefficients modulo reparametrizations, and 
there are 5d + 1 equations, which happens to be the same number. The number of rational 
curves grows exponentially with the degree. There are for example 609250 degree 2 curves and 
704288164978454686113488249750 degree 10 curves. Furthermore, we can consider arbitrary 
D4-charge N. Requiring such a degree N 4-cycle to contain a degree d 2-cycle freezes Nd + 1 
moduli, out of 5Af3 + 25jV . So at large N the 4-cycle can "store" up to degree d rational 
curves by switching on the appropriate fluxes. If d is sufficiently large so that the number 
of degree d curves Nd ~ e Kd is much larger than this number, then this leads, naively at 
least, to up to ~ e"f N2 susy configurations. Switching on one degree d flux quantum increases 
<?o by + 1, and q 2 by ±d depending on the sign of the flux (see section 4.3 of [107] for 

37 Actually we should shift F by the "half-flux" \D X , because for N = 1, S is not spin [111, 112]. This 
ensures proper charge quantization. We will ignore this here. 



73 



a derivation). So at large N the energy in the maximal amount of stored curves is 

9jlU ~ n ~ _X i 5N^_ dN ^ _5iV£ _j_ 5N^_ 57V^ 38 
2?r ~ W ~ 24 ' 6d 2 ~ 24 ~ r 12 24 " 

Although the above discussion is for weak coupling, there is strong evidence that some of 
this structure carries over the strong coupling black hole regime, in the form of the structure 
of multicentered black hole bound states (cf section 5). We refer to [107] for more discussion 
on this. 

4.5 Counting super symmetric states 

Although the above construction is explicit, and the number of ground state configurations 
that can be built in this way is huge, it is still only a small subset, especially at large go- 
In this regime the generic ground state configuration is isolated, and their total number and 
distribution over the moduli space can be computed using the statistical methods to count 
flux vacua developed in [36, 113, 114] and reviewed in section 6 of [7]. Those methods map 
the classical critical point counting problem, after making a continuum approximation, to 
the computation a supersymmetric finite dimensional integral, essentially a finite dimensional 
version of what we did is section 3.8. In the present case we already have a supersymmet- 
ric quantum mechanics to start from. Not surprisingly, it is closely related to the effective 
supersymmetric quantum mechanics used in the statistical approach. 

In what follows we will consider the problem by applying the general machinery of su- 
persymmetric quantum mechanics reviewed is section 3. This clarifies and complements the 
results obtained by the methods reviewed in [7] (and applied to the case at hand in appendix 
G of [107]). In particular it allows us to go beyond the continuum approximation. The com- 
putations will get a little technical, for which I apologize. I chose to include them here because 
it may be useful for some readers to see a nontrivial example worked out, and because as far 
as I know this has not been done in the literature (the treatment has some overlap with [107], 
but we will not use S-duality as an input, but rather derive it directly from the susy quantum 
mechanics) . 

From our general consideration of supersymmetric quantum mechanics in the previous 
section, it is not hard to come up with a formula for the Witten index Q(q) as a weighted 

38 Notice in particular that the energy becomes positive before all moduli are frozen. If what we were 
studying were not D4-branes but D7-branes in IIB orientifolds, this would basically mean that we cannot fix 
all moduli in this way while respecting tadpole cancelation, going against common genericity arguments. This 
is not conclusive though since this construction certainly does not exhaust all possibilities. 
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Figure 13: Artist impression of a smooth abelian D4 configuration with U(l) magnetic fluxes turned on 
(flux density represented by colors). 

sum of Euler characteristics of critical point loci M. s of the superpotential Ws, summed over 
fluxes/momenta S with the given total charge q: 

n(q) = J2(-^ ds x(M s ). (4.73) 

s 

The sign (— l) d ~ ds appears because m transversal dimensions contribute m (the Morse index) 
to the fermion number. This is a useful formula for supersymmetric configurations with small 
DO-charge, but not for large DO-charge. We will try to extract the large charge asymptotics 
from the partition function. 

The supersymmetric partition function generating the Witten indices £l(q) in each con- 
served charge sector is, with H as in (4.61)-(4.63): 

Z((3, g„ <P) = Tr (-if e ^ H ^ = ]T Q(q) e "^(«>+( j A ^) a )-^(*°*+^^) . (4.74) 

<? 

Introducing rescaled RR potentials 

(p = P<f>, (4.75) 
we obtain the following expression for the Witten index £l(q): 

tt(q) = lim / dtp e 27Up - q Z(/3, <p) , (4.76) 
p-*oJ 

where the (p integrals in principle are over a unit imaginary interval, e.g. from — i/2 to i/2. For 
computational purposes it is more convenient to work with integrals over the entire imaginary 
axis, because it allows easy implementation of contour integration techniques. Extending 
the integration ranges in the integrals above only produces trivial additional delta-function 
factors, which in the end can be stripped off manually. 
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We will now try to compute the Witten indices from the Euclidean path integral repre- 
sentation of Z in the limit (3 — » 0, analogous to the derivation leading to (3.36): 

Z = J VzVi)V$e~ s[z ^ A , (4.77) 

with periodic boundary conditions on all fields: z a (r + (3) = z a (r), ip a (r + (3) = ip a (T), 
$*(t + (3) = i? l (r) + 27m*, with ri 1 e Z. By Wick rotating (4.66) to rescaled Euclidean time 

t = -i/3r, (4.78) 

we obtain the Euclidean action 

S= [ dT(^\z\ 2 + ^D t ^ + ^R^^-2nu^ + —(^ + 27Tie)l-—(^ + 2m£) 2 \ 
Jo \9sf3 7T 24 4nv v /+ 4nu v /_ / 

(4.79) 

where we suppressed indices for clarity, denoted £ = L//3, i.e. 

t = (p* - ifi i^VdlT + i/3 i^VdlT , (4.80) 
and introduced "light cone" coordinates 

M =A + ^°, ^^A-^ . (4.81) 

As in (3.35), in the limit — > 0, the path integral localizes to constant paths (z,ip), with the 
fluctuation determinants for 2 and ip canceling each other. The path integral over $ does not 
localize to constant paths, since the $ kinetic terms remain finite when /? — ?• (for (p° ^ 0). 
The d path integral is essentially the partition function of a free particle on a 62-dimensional 
torus. Making the sum over different winding sectors $ l (r + 1) = -$ 4 (r) + 27m* explicit, this is 

J2e-^ n+ie) + + ^ n+ie) - [v#e-^ti dT W+-^-K (4.82) 

n J 

The first factor corresponds to the action weights for the classical straight line trajectories in 
each sector, and the remaining path integral corresponds to the fluctuations from the classical 
paths and thus has strictly periodic boundary conditions d l (r + 1) = $*(r). This is just the 
standard Euclidean free particle propagator integrated over the $-torus: 

| We ^/;^M-^ = u- h */ 2 v- b H\ (4.83) 

By making use of (4.52) and the various orthogonality and (anti-) self-duality properties dis- 
cussed there, as well as the projectors introduced in (4.57) to write explicit expressions for 
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n + and n_ = n — n + , the sum in (4.82) together with the 4-fermion curvature term in (4.79) 
reduces after some tender and care to 

n 

where we separated out a factor in the summand that is independent of the dynamical variables 
z, i/j, in which we also introduced the notation Xj = J l Xi for the projection of X in the 
direction of J, and Xj_j for its orthogonal complement (w.r.t. 77^). Now consider the z,ip 
integral of the other factor: 

d 2d z rfV d d $ e-^> n '\ 2 -^ C^n"+cc)-i« r^) (4g5) 

This integral has the typical zero dimensional supersymmetric form with superpotential W n (z) = 
nTL^z). It localizes on the critical point locus Ai n of W n , meaning the Gaussian approxi- 
mation to the integral is exact. 39 The Gaussian integral for quadratic fluctuations normal to 
A4 n produces a factor for each complex normal direction, while the corresponding normal 
fermions produce a factor Together this gives a factor — for each normal direction. The 
tangential directions to M. n remain integrated over, and the corresponding tangent fermion 
integral gives for each complex tangent dimension a factor —^^-R. Finally, there will be 
an overall path integral normalization factor (2ng s (3)~ d as in (3.35). Altogether the above 
integral thus reduces to (— l) dn e{M. n ) v d u~ d , where d n = dim.M„ and e(A4 n ) is the Euler 
density on Ai n , as in (3.36). The sign is physically meaningful; we expand on this below. 
Putting everything together, we find the generating function: 

Z = {u/v) 1 ' 2 u- b *' 2 e 2 ™% J2(-V dn Xn e -^ n+i ^ + ^ n+i ^ , (4.86) 

n 

where we defined the differential geometric Euler characteristics 

Xn= [ e(R), (4.87) 



and we used d = b2+ 2 ~ 1 ■ Note that since by definition Q n is nonzero only if (n±j) + can be 
made to vanish, all terms in this series have n 2 ±J < 0. 

39 Actually, because n can transform by a monodromy when going around a loop in moduli space, this only 
holds up to boundary terms for a given n, but these cancel between different values of n. To see localization 
more directly, note that since the Witten index is invariant under a change of g s , we can make the integrand 
as sharply peaked on M. n as we want, with vanishing boundary terms when g s — > 0. 
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Defining the dual charges qA = DAin\ q = — £ — ^n 2 , we can also express this as 



Z = (^)V2^ 2 /2 e 2^(n-l)-^i + ^^ ( _ 1) d^ fi( ~ )e -^ g -o--(^i)e + ^^^--¥^^ , (4.88) 

Q 

with as in (4.73) Q(q) = 'Yl, n {—^) dn d Xni summing over the values of n with the given charges 
q. Comparing this to the original (4.74), which can be written as 

Z = e~ 2 ^ uq °-^ u+v)q +- 2 ^ , (4.89) 

we see that what we have shown is roughly that Z is a modular (Jacobi) form under the 
modular transformation u — > 1/u, v — > 1/v, transforming similar to an ordinary theta function 
(although what we have here is much more nontrivial than a theta function, due to the 
nontrivial moduli dynamics). The duality exchanges winding and momentum modes, and 
there are many ways of understanding it: electromagnetic duality of the D-brane theory, S- 
duality of the parent string theory, modularity of the parent CFT in M-theory, T-duality of 
the torus, a version of Poisson resummation, etc. The modular transformation relates the 
large and small flux / DO-charge regimes. 

In any case, we can now extract the large go asymptotics of Q(q) from Z. For simplicity 
we will consider here the case qA = 0, but this is easily extended. Performing the Gaussian 
integrals over the ip A in (4.76) is easy. We get 



/too 
dcp 1 --- d<p B * Z(u, v, if) = cu- k e 2 ™^ 
■ioo 



Xne^. (4.90) 



Here k = ^zJh. w jth 6 2 — dimiJ^S), B 2 = dimif 2 (X), and n± is the component of n 
orthogonal to all cp = ^ a Da-, i-e. the component of n orthogonal to the pullback of H 2 (X) 
in ff 2 (£). Explicitly n\ = n l — (n j D A )D A , where D A = D AB D B , D AB being the inverse of 
Dab = DabcN , with Dabc defined in (4.16). The constant c is the Gaussian determinant 
factor: c = (det Dab)~ 1 ^ 2 - For the example of the quintic we have k = (5iV 3 + 50iV — 3)/2, 
Dn = 5N, c = ^fi- The sum over n is trivially divergent because adding an arbitrary 
m = 7ji a Da € H 2 (X) to n does not alter n±. From the point of view of the integral over 
ip A this is due to the redundancy (p A — > (p A — im A , n — > n + m A DA, which we introduced 
ourselves a little earlier when we extended the integration domain from an interval to the full 
imaginary axis. The upshot is that we can take this into account simply by restricting the 
sum over n G if 2 (S,Z) to a sum over the quotient n G H 2 (Y,, Z)/iJ 2 (E, Z). 
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The go ~~ > 00 asymptotics of Q will be captured by the u — > behavior of Z. In this regime 
all n 7^ exponential corrections can be dropped. The (p° integral can then be evaluated 
exactly by closing the contour. There is a pole of order k at u = i.e. (p° = —/3 < 0. If 
q' Q = g + |j < we must close the contour to the right and the integral vanishes. If q' > 0, 
the contour is closed to the left, and the residue theorem gives 

fifa, 0) - 2n\c\ (-l) d X (M) ^0^- ~ (^) * . (4.91) 

The last approximation is valid for A; 3> 1. It is obtained as a Stirling approximation, or 
equivalently as the saddle point approximation of the integral. The saddle point lies at ip° = 
27iq' /k. This implies we need at least q' 3> k to justify dropping the n 7^ terms, or 
equivalently q ^> \ ~ 

This expresses the number of supersymmetric flux configurations as the volume of a 2/c- 
dimensional shell of radius squared 2q' , which loosely (but not literally because 77^ is indefinite) 
can be thought of as the shell in flux space for which the DO-charge — |j — ^rf^SiSj equals 
go and the D2-charge is constrained to be zero. This reproduces the large go asymptotics 
for the number of flux vacua found in various contexts in [36, 113, 114, 7, 107], based on an 
approximation in which fluxes were taken to be continuous. 

Note that xi-M.) = J M e(R) is the Witten index for the pure D4 without the flux degrees 
of freedom. Since this space is topologically CP d , one may expect x(-M) = x(CP d ) = d + 1 ~ 
N 3 /6. This is not obviously correct because of possible singularities in Ai, but arguments 
were given in [107] that this is nevertheless the correct identification. The sign factor (— l) d 
has a physical meaning: In the large go limit, almost all supersymmetric configurations will 
be isolated critical points of the superpotential Ws(z). Isolated critical points lead to susy 
states with fermion number equal to their Morse index, which here always equals d because 
the superpotential is holomorphic. This explains the (— l) d . The sign (— l) dn in the terms 
with n/0 further suggests (but does not prove) that these correspond to the contributions 
from non-isolated configurations, with d — d n residual moduli. 

It is of course easy to get exponentially large numbers of flux configurations out of these 
estimates, even for modest charges. For example for the quintic with, say, iV = 5, go = 
365, we get Q ~ 10 500 , widely considered to be a very large number. (Similar estimates in 
the context of mathematically very similar constructions of type IIB string vacua with D7- 
branes form the basis for suspicions that there exists a staggeringly huge landscape of string 
vacua sweeping out a for all practical purposes dense set in parameter space, giving a simple 
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solution to the cosmological constant problem but obliterating hopes of top-down predictivity 
[31, 32, 33, 34, 36, 113, 114, 39, 7).) 

On the other hand note that when 2Tceq' < k, the above estimates give exponentially 
small estimates. Clearly then, in this regime, it must be that the terms dominate the 
degeneracy. If it is indeed true, as suggested above, that contributions can be identified 

with non-isolated configurations, this means that in this regime, non-isolated configurations 
become entropically dominant. 

For q' <C k we could try to build up the spectrum along the lines of section 4.4, and this 
could in turn be used to infer the corrections in the regime q' S> k. However in the intermediate 
regime q' ~ k, the system is extremely complex. 40 Trying to count degeneracies in this regime 
may be like trying to compute the boiling point of water from first principles. In the case 
of water we can just measure the boiling point and be done with it. We let nature do the 
computation for us. In the case of D-branes, this would not appear to be an option. However, 
there is something analogous, provided we enlarge our task to the problem of counting all D4- 
D2-D0 bound states, not just those corresponding to abelian flux configurations of the D4. 
In that case, we can construct the corresponding BPS black hole solutions, and simply read 
off their Bekenstein-Hawking entropy S = \ogQ(q). We let gravity do the computation for 
us! To include all entropically relevant D4-D2-D0 bound states we need in particular consider 
bound states with localized DO-branes. We turn to this next. 

4.6 Bound states with DO-branes 

So far we have only considered D4-brane and U(l) worldvolume flux degrees of freedom. A 
simple but entropically important extension is to include bound states with pointlike DO- 
branes. For a D4 wrapped on a fixed smooth E, the number of such bound states is easy to 
compute [115]. We can think of the DO-branes as a gas of noninteracting particles. DO-branes 
form bound states among themselves of arbitrary DO-charge k (the easiest way to see this is 
to consider their uplift to M-theory, where they are KK modes along the M-theory circle, k 
being the KK momentum), so each particle in the D4-D0 gas is characterized by a charge 
quantum number k > 0, as well as by the supersymmetric 1-particle state \pa) it occupies. 

40 This happens to be the relevant regime for applications to type IIB string theory vacua [7, 105]. In this 
context space-localized D4-branes are replaced with space-filling D7-branes, carrying U(l) fluxes inducing 
D5- and D3-charges. The D3-tadpole cancelation condition dictates that the total D3-charge vanishes, hence 
9o = 24 ~ an d 27r fc <?0 ~ ~ 1.42, barely above the threshold. 
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Figure 14: Artist impression of a smooth abelian D4 configuration with U(l) magnetic fluxes and mobile 
DO-branes bound to it. 

As we have seen in section 3.5, these 1-particle states are given by the harmonic differential 
forms on E. The label p = 0,1,2,3,4 in \pa) indicates the fermion number (form degree) 
and a = l,...,b p . By the usual Fock space construction, an arbitrary multi-particle state 
can be represented by specifying occupation numbers nk, p , a , which take values in Z + if p is 
even and in {0, 1} is p is odd. In the absence of fluxes, the total DO-charge of the D4-D0 
state is Q = — ^ + Ylkpa ^ n kpa and the fermion number is F = YlkpaP n kpa- Here we are 
still ignoring the moduli degrees of freedom of E, considering those frozen for the time being. 
The generating function for the degeneracies cIqf of supersymmetric ground states of total 
DO-charge Q and fermion number F is then 

^^dQpy F t® = t~24 yT.kpc,P n kpa -j-T.kpa kn kpc 

QF {n ka } 

= r& YlJ2y np t nk 

kpa n 

= ft tt IWI + »*>'■ (4 92) 

A generating function for the Witten indices ^do(Q) = J2f(~ ^) F ^qf is obtained by setting 
y = — 1 in the above: 

]T n m (Q) e = r^ - t k y* = nit)-*, (4.93) 

Q k 

where % = X! P ( — l) p ^p i s the Euler characteristic of E, and r](t) is the Dedekind eta-function. 
If we ignore the fact that there can be nontrivial interplay between the D4 and DO moduli, 41 



41 Ignoring this is actually not a good idea at small go, but at large qo it is a good approximation since most 
D4 flux configurations will be isolated anyway in this regime. 
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adding the DO-branes to the D4 simply has the effect of multiplying the partition function 
Z in (4.76) by e~ 27TU 2iri(e~ 27ru )~ x . Using the modularity property of the eta function, we 
furthermore have 

v {e- 2 ™)- x = u x/2 r](e-^ Y x . (4.94) 

Combining this with the D4 partition function (4.86), using \ = b 2 + 2, and for notational 
simplicity at the cost of slight formality setting v = —u = — <£>° in (4.86), we get, within our 
no D4-D0 interplay approximation: 

n 

2 9 
2-7T i x \ i — r 2nk v , 2tt r n , ■ \ 

= ^^i^--) _ e -^r)-x J2(-l) dn x n e^ { ~ + v) . (4.96) 

k n 

Again we see that Z transforms as a modular form under ip° — > Notice in particular 

that the large positive and negative weights of the D4 resp. the DO-part have canceled. 

To extract the large go asymptotics, we retaining only the leading terms at small positive 
<£>° . After integrating out (p: 

Q(q) « J V # +1 e (-l) d X (M) . (4.97) 

where q = q + — — ^ Ms -. The saddle point of this integral is at (p° = J with value 

n(q)^e 2n ^r. (4.98) 

We recall that \ = D ABC N A N B N C + c 2 , A N A , so this is an explicit formula for the index 
function of the charges. The logarithm of this expression can therefore directly be compared 
to the corresponding black hole entropy. 

This is reliable when the saddle point value of ip° is small, i.e. go 3> X- When the latter is 
not the case, the saddle point value of tp° is not small, and so, as in the case without DO-branes, 
even within the model limitations we have made (such as ignoring singular or nonabelian D- 
branes), the above approximation for Q(q) is not necessarily reliable. We will see in the next 
section that the agreement with the black hole entropy is excellent when g ^> x- I n f ac t 
for the leading order matching in this regime, we do not even need to include the fluxes and 
moduli degrees of freedom, and a much simpler derivation is possible [93, 109, 115]. The 
match does get better at subleading order when the D4 degrees of freedom are included. But 
when g < X-, the correct entropy is not anywhere near (4.98). In this regime, the flux and 
D4-moduli degrees of freedom dominate the entropy. 
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4.7 Some extensions 



With the goal of introducing in detail a complex, glassy system that occurs naturally in string 
theory and has a precise geometric description and well-controlled holographic counterparts, 
we have given a fully explicit construction of the supersymmetric quantum mechanics de- 
scribing the ground state sector of wrapped D4-branes bound to DO-branes, and discussed 
explicit constructions of supersymmetric ground states and their counting in some detail. 
Some of this was a review of parts of [107] and the ideas and results used there, with some 
simplifications and some points worked out more explicitly. Our discussion was necessarily 
incomplete in scope. Rather than try to give an overview of the huge related literature, let 
me mention a few immediate extensions that could have directly followed this part, with a 
small, non-representative, biased sample of possible starting points for further reading: 

1. The lift to M-theory as an M5 wrapped on S 1 x S, reducing to a (0,4) CFT. This is in 
some ways a more natural framework as it unifies D4 and DO degrees of freedom and 
allows for a holographic dual description beyond the ground state sector [109, 116, 117]. 

2. Quiver quantum mechanics: these are Id gauged linear sigma models, providing a simple 
but very rich class of models describing complex D-brane systems. The field content is 
represented by quiver diagrams, nodes being partonic branes and arrows open strings 
between them [118]. Quantization allows explicit interpolation between the weakly cou- 
pled geometric regime and the strongly coupled black hole bound state regimes [119] 
(see also below). 

3. A more refined description of D4-D0 bound states that takes into account the interplay 
between D4 and DO moduli, nonabelian degrees of freedom, and so on. One conceptually 
simple but efficient approach is the brane-anti-brane tachyon condensation picture [120, 
121, 122, 123, 124, 107, 125]. 

4. Counting of ground states away from the regime go ^ X- Relation to Gromov-Witten, 
Gopakumar-Vafa and Donaldson-Thomas invariants and the OSV conjecture [104, 126, 
127, 128, 129, 130, 107, 131]. 

5. Wall crossing phenomena: e.g. the possibility for a D4 to split into a D6 and an anti-D6 
away from the strict large volume limit [132, 121, 133, 119, 107, 134, 135]. 

6. Applications to global aspects of IIB model building [125, 136] 
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It would certainly be useful to have a fully explicit example in which in particular the period 
vectors can be computed exactly. This may be possible along the lines of [137]. 

We end with some references to key papers in the history of the subject. The idea that 
branes can be wrapped on nontrivial compact cycles to obtain charged particles in four dimen- 
sions appeared first in [138] (and shown to be necessary for the consistency of string theory in 
[139]). The discovery of the perturbative string description of D-branes [140] made it possible 
to quantize and count states at weak coupling, leading to the computation of the Bekenstein- 
Hawking entropy in [93]. Many of the ideas in the more systematic development of the physics 
and mathematics of wrapped D-branes in Calabi-Yau manifolds, including stability and wall 
crossing, originated in [121, 141], with precursors in [142, 132]. 

5 Black hole duals 

5.1 Introduction 

In the appropriate regime, D-branes wrapped on compact cycles manifest themselves as black 
holes in the low energy effective field theory. It turns out that an intricate zoo of stationary 
supersymmetric bound states of such black holes exists, like giant molecules, all of which 
collapse to a D-brane localized at a single point in space in the limit g s — > 0. 

In this final section we summarize the explicit solutions of such black hole bound states, 
and briefly discuss their entropy compared to the microscopic picture developed in the previous 
section. The complexity of the black hole solutions is correlated with the complexity of the 
D-brane landscapes that arise at weak coupling. We will focus again on type IIA Calabi-Yau 
compactifications. 

5.2 BPS solutions of J\f = 2 supergravity 

The four dimensional effective theory of type IIA string theory compactified on a Calabi-Yau 
manifold A is M = 2 supergravity coupled to B 2 = 62(A) abelian M = 2 vector multiplets. 
There are B 2 + 1 gauge fields A A , A = 0, 1, ■ ■ ■ , B 2 , obtained as in (4.13) by KK reduction of 
RR potentials: C^ 1 ' = A , = A a ADa- The B 2 vector multiplets furthermore each contain 
a complex scalar, obtained by KK reduction of the complexified Kahler form: B + iJ = t A DA- 
They also contain spin 1/2 fermions, but we don't need those. Finally, the theory contains 
massless hypermultiplets too, but they do not affect the solutions of interest to us and can be 
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consistently put to constant values. We will work in this section in units with the 4d Newton 
constant Gjy = 1. 

Wrapped D-branes manifest themselves as charged point particles in four dimensions, 
sourcing the vector multiplet fields. We have both electric (D2/D0) and magnetic (D4/D6) 
monopole charges. In general, the lattice L of magnetic-electric charges F carries a fundamen- 
tal symplectic product, which in a symplectic basis has the canonical form 

(r,f) =r A f A -r A f A . (5.1) 

Upper indices denote magnetic and lower indices electric components. In the IIA case at 
hand, the electric charges are the DO and D2 charges Qq and Qa, and the dual magnetic 
charges are the D6 and D4 charges N° and N A . A2 In terms of these, the symplectic product is 
((N, Q), (N, Q)) = N°Q + N A Q A - Q N° - Q A N A . Integrality of this product is equivalent 
to Dirac quantization. 

In the weak string coupling regime (4.48), the wrapped D-branes are well described as point 
particles moving in flat M 3 , interacting with each other primarily through the lightest stretched 
open string modes. The coupling g 4 = g s /VVx oc £ p /£ s is in a hypermultiplet and can 
therefore be tuned at will. When it gets larger, excited open string modes become important, 
until eventually the interactions are better described by massless closed string exchange, i.e. 
graviton, photon and scalar exchange. This is the regime in which the supergravity description 
becomes valid. When the charges are large, the wrapped D-brane states manifest themselves 
as weakly curved black hole solutions. 

A single centered BPS (i.e. supersymmetric) solution to the equations of motion is neces- 
sarily static and spherically symmetric, with a metric of the form ds 2 = —e 2U dt 2 + e~ 2U dx 2 , 
and all fields functions of r = \x\ only. The BPS equations of motion take the first order 
gradient flow form typical for supersymmetric solutions [143, 144]: 

U=-e u \Z\, (5.2) 
i A = -2e u g AB d s \Z\, (5.3) 

where g A g is the metric on the vector multiplet moduli space, the dot denotes derivation with 
respect to r = 1/r, and Z(T,t) is the central charge of the magnetic-electric charge T in a 

42 Usually these are denoted by P° and P , but to be consistent with the previous sections we use iV and 
N A instead. 
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where 



background with moduli t A . It is a complex function on the vector multiplet moduli space, 
holomorphic up to a normalization factor, and linear in the charge vector V: 

V A = X\ V A = -^, X A = t A X\ A = 1,...,B 2 . (5.5) 

Here F(X) is the prepotential of the M = 2 theory, which determines all couplings and metrics 
in the 4d action. In general it is a locally defined holomorphic function, homogeneous of degree 
2. For our type IIA theory, it takes the form 

F(X) = -\d abc X a X b X c + ■■■ (5.6) 
6 

where D AB c was defined in (4.16) and the ellipsis denote string worldsheet instanton cor- 
rections, which are exponentially suppressed in the Kahler moduli J A = lmt A and therefore 
negligible when the Calabi-Yau is large. Dropping those, (5.4) boils down to 

1 — ( l -N«D ABC t A t B t c - l -N A D ABC t B t c + Q A t A + Q ) . (5.7) 

8D ABC J A J B J C \ ' V 6 2 / 



The central charge gets its name because it appears as a charge commuting with everything 
in the M = 2 supersymmetry algebra. Its absolute value equals the lowest mass a particle of 
charge T can possibly have in a background specified by t. This bound is saturated for BPS 
states. Its phase determines the supercharges preserved by the BPS state. Two BPS objects 
are mutually supersymmetric if their central charge phases line up. 

The gradient flow equations (5.2) drive \Z(T, t) \ to its minimal value |.Z*|. 43 If this minimal 
value is zero at a nonsingular point in moduli space, no solution exists. If it is nonzero, we 
get a black hole with near horizon solution t = const. = £*, e~ 2U ^ = as can be checked 
directly from the above BPS flow equations. This describes AdS2 x S 2 with S 2 horizon area 
A = At[\ZJ\ 2 and therefore Bekenstein-Hawking entropy 

5(r) = 7r|^(r)| 2 . (5.8) 



43 Although all regular critical points of \Z\ are isolated local minima [145], in the presence of singularities 
at finite distance in the moduli space, there can be multiple basins of attraction, but we will ignore this here, 
as it does not occur in the large volume approximation. 
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Figure 15: Warp factor e u for some three centered D6 — D6 — D4 bound state, for a choice of charges that 
doesn't really matter because this figure is primarily ornamental. 

When the D6-charge N° is zero, the minimization of \Z\ is straightforward and there is a 
simple closed form expression for S(T) (still in the large volume approximation) [146]: 



where Dab — DabcN and D is its inverse. Note that to leading order in the large charge 
limit, this is exactly the same as the microscopic result (4.97)! The subleading term can be 
reproduced macroscopically as well, from the Wald entropy in the presence of an R 2 term 
[147]. When N° is not zero there is no general closed form solution for S(T), except when Qa 
is chosen to be proportional to DabcN b N c . In that case we define for some fixed K A > 
the charge parameters (n°,n,q,q ) by (N°, N A , Qa, Qq) = (n° ,n K A ,q (K 2 )A,qo K 3 ), where 
K 3 is a short for D A bcK a K b K c and (K 2 ) A for D A bcK b K c . Then we can write 



As we will see below, once the entropy function S(T) is known, it is easy to write down fully 
explicit expressions for all fields at all points in space. 

Quite remarkably, these theories also have multicentered, supersymmetric, stationary black 
hole bound states, like giant molecules. These are genuine bound states in the sense that 




(5.9) 



S = - K 3 J 3 n 2 q 2 - 8 n° q 3 + 6 n 3 q° - 18 n n° q q° - 9 n° 2 q° 2 . 



(5.10) 
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the centers are constrained by nontrivial potentials, generated by scalar, electomagnetic and 
gravitational forces [133, 148, 149, 150, 117, 107]. They have a metric of the form 



ds 2 = -e 2U (dt + Uidx'f + e~ 2U dx 2 , (5.11) 

where U and u depend on x, and they are fully characterized by harmonic functions if A , H\ 
with sources at the positions Xi of the charges Tf. 

H = T,]^-\+ h > ( 5 - 12 ) 

i 

where the constant term h is determined by the total charge T and the asymptotic moduli: 

h = -2Im(e-"V)| r=00 , a = argZ(r), V = V/\(V, V)\ l/2 , (5.13) 
and V is as in (5.5). 

As mentioned earlier, once the entropy function S(T) is known on charge space, the com- 
plete solution is known [150], simply by substituting H(x) for r in S and its derivatives: 



e-™ = \s { H), ^=i^||^(* + W ) + ^ on , ^ = |^|#. (5.14) 

The one-form v4^ on is the vector potential for a system of Dirac magnetic monopoles of charge 
located at the positions Xj. The off-diagonal components u> of the metric are the solutions 

to 

V x u = (VH, H) , (5.15) 

where V is the flat space gradient. This equation implies an important integrability condition: 
V • (V x to) = =>- (V 2 H,H) = 0, from which, using V 2 |4| = — 47r<5 3 (x), we get for every 
center i a condition: 

Eifrl = "< r ^> = 2Im ( e ~ ia ^ r ^))Loc • ( 5 - 16 ) 
j j 

This imposes constraints on the positions xf, solutions are the BPS equilibrium positions of 
the black holes in their mutual force fields. For the 2-centered case we have in particular 

| il -4| = MM| | (5 . 17) 



J 2\ 

where = Z(Ti,t)\ r=00 . When the right hand side goes from positive to negative through a 
wall where Z\ and Z 2 line up, the equilibrium separation diverges and the bound state decays. 
This is the supergravity incarnation of the wall crossing phenomenon. 
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Figure 16: (from [107]) Attractor flow tree in the i-plane for a particular bound state of 7 flux-carrying 
D6-branes and 7 flux-carrying anti-D6-branes. The flow starts at the yellow dot at the top and goes down. 
The iterated structure can be thought of as representing clutering of individual centers. In the case at hand all 
attractor points lie on the boundary of moduli space, which is due the special choice of charges. The symmetric 
form is also special to our particular choice. The main motivation to pick those charges was esthetics and 
dramatic effect. More generic examples can be found in [107]. 

For completeness we also give the solution to (5.15) [150]. For a 2-centered bound state 
with centers at (0,0, L) and (0,0, — L), up to residual diffeomorphism gauge transformations 
t ->■ t + f(x), u ->■ uj — df: 

(T T ) { L 2 — r 2 \ 

1^2/ ^^ + 1 _ COS 1+CO s0 2)d6. (5.18) 

2L \(L 4 + r 4 -2L 2 r 2 cos2#) 1 /2 2 J Y K J 

Here (r,<p,8) are standard spherical coordinates centered at the origin and d±, 6 2 are angles 
of spherical coordinates centered at the two particle positions. The integrability condition is 
equivalent to the absence of physical singularities along the z-axis. Because (5.15) is linear, 
the solution for more centers is obtained by superposition. 

The solutions are generically stationary but not static: they have intrinsic angular mo- 
mentum, given by 

J = -^2(Ti,Tj) Xij , Xij = — J ' . (5.19) 

^ . \Xi Xj 

Kj 

This generalizes the well-known electromagnetic field angular momentum sourced by monopole- 
electron pairs (generated by the E x B Poynting vector circling around the electron-monopole 
axis) . 
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The integrability conditions are a necessary but not sufficient condition for the existence of 
solutions. One also needs to check if the entropy function S(H) is positive everywhere. This 
is in general difficult to do. In [133, 107] the conjecture was formulated that the existence 
of solutions is equivalent to the existence of the (much easier to analyze) attractor flow trees. 
An attractor flow tree is built out of single attractor flows. The tree starts at the background 
value of the moduli and terminates at the attractor points of the constituent charges. Each 
edge E of an attractor flow tree is given by a single charge attractor flow for some charge r^. 
Charge and energy are conserved at the vertices, i.e. for each vertex splitting E — > (Ei,E 2 ), 
T E = T El + Y E2 and |Z(r^)| = |Z(r Sl )| + |Z(r^ 2 )|. The last condition is equivalent to 
requiring the vertices to lie on a line of marginal stability: argZ(r^ 1 ) = axg Z{Ye 2 )- Flow 
trees can be thought of as "skeletons" for the supergravity solutions, or as giving a recipe 
for adiabatically assembling or disassembling multicentered bound states. This gives a strong 
physical argument for the conjecture. We refer to section 3.2.2 of [107] for a more detailed 
discussion. 

5.3 Entropy 

Comparing the macroscopic and microscopic D4-D0 entropies (5.9) and (4.97), we find perfect 
agreement to leading order when Q ^> iV 3 . The sub leading correction can be reproduced 
from the Wald entropy in the presence of an R 2 correction term in the 4d effective gravity 
theory [147]. In fact, in the same regime, an infinite series of corrections matches [104, 107, 
128, 130, 129] 

However, as was clear from the discussion in section 4, microcopically things get much 
more complicated when Q becomes of the order of x ~ the D4-degrees of freedom start 
getting important and eventually dominate the entropy. On the black hole side things get 
much more complicated too in exactly the same regime: Multicentered solutions appear, and 
at some point they start dominating the entropy [107]. 

In particular one clearly exits the easy regime where Qo ^ A^ 3 whenever all charges 
are uniformly scaled up sufficiently. Thus, more precisely, when charges T are scaled up as 
T — > AT, there may exist for example two-centered solutions with horizon entropy scaling as 
A 3 , while the single centered entropy scales as A 2 . 

Let us give a concrete example of this phenomenon, taken from [151]. Consider the Calabi- 
Yau X = Tl x T| x T 3 2 , a product of three two-tori. Let £ be the 4-cycle (Tf x T|) + (Tf x 
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T|) + (Tg x T 2 ) and let a be the 2-cycle T 2 + T 2 2 + T 3 2 . Then the entropy function of a charge T 
corresponding to wrapping N° D6-branes on X, N D4-branes on £, Q D2-branes on a and Qo 
DO-branes is given by S(T) = n^AN 3 Q Q + 3N 2 Q 2 + QN°NQQ - AN°Q 3 - (N°Q ) 2 . Now 
consider a total charge 

T = (N°, N, Q, Qo) = A(0, 6, 0, -12), (5.20) 

in a background in which the area of each T 2 equals v. Then for any v, there exists a single 
centered solution with horizon entropy 

5i = 72v / 2ttA 2 . (5.21) 

However, there also exists D6-anti-D6 two-centered black hole bound states, for instance with 
charges 

T 1 = (1, 3A, 6A 2 , -6A), T 2 = (-1, 3A, -6A 2 , -6A). (5.22) 
provided v > \/l8A. The constant terms in the harmonic functions (5.12) are h = (0, h=, 0, — 



and the equilibrium separation (5.17) at large A is \xi — x%\ ~ T^p^t • ^ ne resu lt m § metric 
is well defined, with S(H) real and positive everywhere. The two centers have equal horizon 
entropy, summing up to a total entropy 

S 2 ^l2nVSA 3 , (5.23) 

which is indeed parametrically larger than the single centered entropy (5.21). 

When v is kept fixed while sending A — > oo, the equilibrium separation diverges and then 
formally goes negative, meaning these 2-centered solutions disappear from the spectrum. But 
the low energy effective quantum mechanics of section 4.3 is supposed to become increasingly 
reliable in the limit of large volume v — > oo (and weak string coupling), so its Witten index 
for given charges must be the total index of all black hole configurations in the limit v — > oo. 
Thus we conclude that at large A, the usual weakly coupled, weakly curved wrapped D4- 
branes are not computing the entropy of a single D4-D2-D0 black hole, but predominantly 
that of multicentered configurations. The single center D4 black hole entropy can only be 
microscopically reproduced from a weak coupling computation in the regime Qo ^ N 3 . 

This occurs in particular also in the lift to M-theory [117], where a precise holographic 
correspondence exists between a 1+1 dimensional (0,4) CFT [109] and quantum gravity in 
AdS 3 x S 2 . It was pointed out there that the solution with the highest entropy in the A — > oo 
limit appears to be a configuration in which all entropy carrying charges has been concentrated 
in one of the D6-charged black holes. 
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5.4 Landscape structure 

In the minimal charge case, namely Qo = — 24 ~ — ^T> ^ ne supergravity solution is a pure 
D6-anti-D6 bound state, where the D6 branes carry flux but no entropic degrees of freedom 
[107]. Quantization of the centers leads to a spin j = d/2 multiplet of supersymmetric ground 
states with degeneracy 2j + 1 = d + 1, in agreement with the microscopic picture where this 
arises as the Euler characteristic of the moduli space M = CF d . The lifted M-theory solution 
becomes pure "spinning" AdS3 x S 2 in this case. 

Some identifications between the microscopic structure of the D-brane landscape and the 
macroscopic black hole solutions have been made (where it should be kept in mind of course 
that the map from D-brane microstates to black hole macroscopic states must be many-to-one, 
due to the existence of black hole entropy). For example switching on fluxes represented by 
holomorphic 2-cycles as in section 4.4 corresponds to "dressing up" the D6 or anti-D6 by halos 
of D2-D0 particles [107], or equivalently in M-theory by adding M2-branes to the north or 
south poles of the S 2 [130, 129]. But the full map at higher DO-charge is not well understood. 

Dynamics near zero energy should proceed through tunneling of charged particles in and 
out of the black hole centers in this picture. This is under study in [3]. 
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